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Preface, 


The  purpose  of  this  study  was  to  examine  the  effects  of 
atmospheric  drag  on  the  long  term  orbit  of  artificial 
satellites  about  the  planet  Mars  using  a  dynamic  atmosphere 
model  in  place  of  a  static  atmosphere  model. 


Existing  literature  on  the  Martian  atmosphere  and 


orbital  equation 

models 

were  examined  and 

then  incorporated 

into  theoretical 

and 

computational  models.  However, 

the 

resulting  models 

are 

only  the  starting 

point  for 

the 

development  of 

more 

sophistocated 

theoretical 

and 

computational  models. 

This  thesis 

establishes 

the 

fundamental  differences,  if  any,  i.n  satellite  motion  due  to 
the  use  of  a  dynamic  atmosphere  model  instead  of  a  static 
atmosphere  modal. 

I  am  indebted  to  my  thesis  advisor,  Capt  Rodn>;y  Bain, 
for  his  suggestions  and  ready  assistance.  He  has  increased 
my  understanding  of  both  the  mathematical  and  orbital 
theory,  as  well  as  providing  a  'relief  valve'  for  all  the 
frustrations  associated  with  the  completion  of  this  thesis. 
I  would  especially  like  to  thank  Richard  P.  Osedacz  for  his 
help  and  moral  support. 
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Abstract 


This  study  examines  the  effects  of  atmospheric  drag  on 
the  long  term  orbital  motion  of  an  artificial  satellite 
about  the  planet  Mars  using  a  dynamic  Martian  atmosphere 
model  Instead  of  a  static  atmosphere  model. 

Present  studies  of  motion  through  the  Martian 
atmosphere  use  a  non-rotating  (static)  atmosphere  model 
which  has  a  simple  exponential  drop-off  with  altitude. 
However,  Mars,  like  other  planets,  is  known  to  have  a 
turbulent,  complex  atmosphere.  As  a  result,  a  new  dynamic 
model  of  the  Martian  atmosphere  uses  data  from  recent  space 
probes  of  Mars  to  simulate  an  atmosphere  that  is  both 
position  dependent,  through  latitude  and  longitude  effects, 
and  time  dependent,  through  a  variety  of  effects.  These 
include  diurnal  and  seasonal  effects,  annual  motion  effects, 
solar  activity  effects,  and  dust  storm  effects.  Nine 
constituent  gases  are  included  in  the  above  model.  The 
model  determines  atmospheric  mass  density  for  the 
calculation  of  atmospheric  drag.  The  existing  static  model 
is  also  used  to  calculate  atmospheric  drag  to  provide  a 
comparison  in  conjunction  with  other  effects. 


xiii 


In  addition  to  atmospheric  drag,  perturbations  due  to 
the  attraction  of  the  sun,  the  aspherical  gravity  field  of 
Mars,  and  the  effects  of  the  solar  wind  are  also  modeled. 

The  resulting  orbital  data  are  then  compared  to 
determine  the  differences  between  the  two  models,  rhus 
showing  the  effects  on  long  term  motion  of  using  a  dynamic 
model  of  the  Martian  atmosphere. 
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AN  INVESTIGATION  OF  LONG  TERM  ORBITS 
ABOUT  THE  PLANET  MARS  USING  A  DYNAMIC 
ATMOSPHERE  MODEL 

lx.  lEtraduc.fci.Qfl 

Current  studies  of  orbital  notion  about  Mars  use  a 
simple  non-rotating  (static)  exponential  model  for  the 
Martian  atmosphere.  However,  the  Martian  atmosphere  is  far 
from  static;  the  entire  surface  of  Mars  has  been  obscured  by 
dust  storms  (37:206).  Recent  data  indicate  the  mass  density 
of  the  Martian  atmosphere  is  heavily  time  and  position 
dependent  (dynamic)  due  to  a  variety  of  perturbing  effects. 
As  a  result,  no  accurate  prediction  of  atmospheric  mass 
density  at  a  particular  position  and  time  in  the  future  is 
possible  (8:1).  However,  using  a  time  and  position 
dependent  model  of  the  atmosphere  over  a  sufficiently  long 
time,  mass  density  can  be  predicted  with  acceptable  accuracy 
since  long  time  intervals  average  out  short  term 
fluctuations  that  cannot  be  modeled  (9:1). 

Current  literature  on  the  Martian  atmosphere  is 
examined  and  then  a  model  for  satellite  motion  through  a 
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"realistic"  dynamic  Martian  atmosphere  is  formulated.  To 
assess  the  effect  of  using  a  dynamic  atmosphere  model,  the 
resulting  orbital  data  from  the  more  realistic  dynamic 
atmosphere  model  was  compared  to  the  orbital  data  from  the 
current  static  atmosphere  model.  Two  advantages  of  using  a 
more  realistic  model  would  be  a  better  guaranteed  lifetime 
for  the  satellite,  required  for  various  NASA  guidelines 
(9:1),  and  a  more  effective  use  of  the  limited  spacecraft 
mass  and  internal  space  in  the  tradeoff  between  fuel  and 
scientific  payloads. 

The  trajectory  propagation  model  includes  perturbation 
models  for  the  aspherical  Martian  gravity  field,  the 
gravitational  attraction  of  the  sun,  solar  radiation 
pressure,  and  atmospheric  drag  at  altitudes  below  1000  km. 


II...  I  he,  Martian 


In  order  to  develop  the  dynamic  atmosphere  model, 
literature  on  the  Martian  atmosphere  was  surveyed  and  then 
used  to  formulate  the  model.  This  section  briefly  reviews 
the  literature. 

SiailArilics  Qf  Mara  to  Earth 

The  key  to  understanding  the  development  of  models  of 
the  Martian  atmosphere  is  to  examine  the  similarities  of 
Mars  to  Earth  and  the  limitations  of  viewing  Mars  through 
the  Earth’s  atmosphere.  Table  1.1  presents  Earth  and  Mars 
information. 

Several  similarities  exist:  1)  the  rotation  rates  are 
nearly  the  same,  2)  both  have  seasons  although  planetary 
tilt  causes  longer  seasons  on  Mars,  and  3)  both  planets  have 
polar  caps  and  clouds.  On  the  other  hand,  1)  Mars  receives 
less  radiation  from  the  Sun  since  it  is  farther  away,  2)  it 
has  no  oceans,  3)  its  atmosphere  is  only  about  0.01  that  of 
Earth,  and  4)  its  polar  caps  are  not  made  from  water. 
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Table  1.1: 

Earth  and  Mars 

Data 

Planet: 

Earth 

Mars 

Mean  Distance  from  sun: 

1.496E8  km 

227,941,963  km 

Eccentricity; 

0.017 

0.093387 

Inclination  to  ecliptic 

:  0* 

1.86641* 

Sidereal  Period: 

365.256  days 

686.9804  days 

Equatorial  Radius: 

6378.2  km 

3396.9  km 

Polar  Radius: 

6356.2  km 

3376.1  km 

Mean  Surface  Gravity: 

978  cm/s" 

372  cm/s2 

Mean  Escape  Velocity: 

11.2  km/s 

5.022  km/s 

Sidereal  Day: 

23  h  56  m  4  s 

24  h  37  m  22.663  s 

Inclination  of  Equator 

to  Orbit:23.5‘ 

24.936* 

Visual  Albedo: 

0.36 

0.159 

(Various  Sources) 

(Various  Sources) 

Furthermore,  due  to  the  orientation  and  eccentricity  of  the 
orbit,  the  Martian  seasons  are  not  the  same  length.  The 
southern  summer  is  short  and  hot,  while  the  northern  summer 
is  long  and  cool,  thus  strongly  affecting  the  polar  caps  and 
the  atmosphere. 

Brief  History  of  QbafirYAiiQaa/Tticoriea 

A  wide  variety  of  seasonal  and  meteorological  data  have 
been  observed  since  Sir  William  Herschel  first  noted  Mars’ 
thin  atmosphere  .  Until  recently,  most  observations  of  Mars 
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were  taken  from  the  surface  of  rhe  Earth  I2U.  However, 

Earth’s  atmosphere  affects  the  data.  Similar  molecules 
contaminate  data,  resolution  is  limited,  and  signal 
strength  is  weakened.  Although  various  means  have  been  used 
to  reduce  the  amount  of  atmosphere  the  signal  travels 
through,  these  also  affected  the  data,  leaving  one 
solution:  space  probes.  A  stronger  signal  and  no  filtering 
atmosphere  mean  the  space  probes  receive  data  too  weak  to 
reach  Earth,  such  as  ultraviolet,  and  have  greater 

resolution  and  accuracy.  They  overturned  many  theories. 

Early  theories  on  the  Martian  atmosphere  reveal  some 

good  deductions  and  some  bad  assumptions.  For  example, 

predictions  of  large  diurnal  temperature  fluctuations  were 

correct  .  However,  the  use  of  terrestrial  assumptions  and 

atmosphere- f iltered  data  led  to  errors.  For  example,  the 

atmosphere,  thought  to  be  mainly  N*2  as  on  Earth,  is  actually 

composed  almost  entirely  of  CO2  ( 5:1-18  ).  Also,  the  mean 

« 

surface  pressure,  which  polariraetric  data  indicated  to  be 
about  85  mbars  (21:23),  is  now  accepted  as  about  6  mbars 
(15:4321).  Further  differences  will  be  mentioned  in  later 
sections . 

Most  Recent  Information/Beliefs  about  the  Martian  Atmosphere 

Major  seasonal  changes  occur  and  are  closely  associated 
with  changes  in  the  polar  caps.  In  early  spring  the 
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southern  polar  cap  has  its  greatest  size.  \3  spring 
proceeds,  the  cap  sublimates  and  begins  to  recede  at  about 
l  of  latitude  every  five  days.  One- third  of  the  Martian 
vear  MtPr  the  cap  has  receded  to  its  smallest  extent, 
len\  hi?  smal  l  remnant  ice  caps  behind  rJ2 :  11*19 ) .  Th* 
t*eh.n  tor  of  the  northern  cap  is  similar,  but  not  ail  me 
details  are  the  same  due  to  the  difference  in  season  length 
and  variation  in  distance  from  the  sun;  the  planet  is  at  a 
greater  distance  from  the  sun  during  the  northern  summer  and 
so  receives  less  heat.  Both  caps  advance  and  retreat  by 
about  35  latitude  each  season.  Mariners  7  and  9  infrared 
radiometer  and  spectrometer  measurements  showed  that  the 
main  polar  caps  are  <,02  (3:404).  Viking  2  showed  the  north 
polar  remnant  to  be  water  ice,  probably  the  largest 
collection  of  water  on  the  planet.  The  south  polar  remnant 
composition  is  still  in  some  doubt.  The  polar  caps  affect 
atmospheric  density  by  causing  the  carbon  dioxide  to  "snow 
out"  of  the  atmosphere  or  to  sublime,  causing  a  seasonal 
variation  of  about  25-30X  (42:3969).  At  the  cap  edge  there 
exist  strong  winds  driven  by  mass  outflow  and  large 
temperature  gradients  during  sublimation  as  well  as 
traveling  waves  (weather  systems)  generated  equatorward 
(49:298).  As  a  result,  changes  in  polar  temperature  have 
profound  effects  on  atmospheric  dynamics  (5:124).  Oust 
storms  affect  the  polar  caps  by  filtering  the  sunlight  and 
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by  contributing  dust  particles  to  the  polar  deposits, 
lowering  their  albedo  (38:270). 

Oust  is  always  being  injected  into  the  Martian 
atmosphere  (38:259).  However,  it  is  greatest  at  times  close 
to  oroltal  perihelion  and  summer  solstice  in  the  southern 
hemisphere.  There  are  about  iOO  local  dust  storms  each 
year,  mostly  during  the  southern  hemisphere  spring  (38:262). 
Once  or  twice  a  year  these  develop  into  global  dust  storms, 
usually  starting  in  certain  preferred  locations.  These 
storms  are  never  quite  the  same,  so  a  long  term  average  must 
be  used  to  quantify  their  effects  (49:294).  Viking  showed 
that  measured  amplitudes  of  thermally  driven,  planetary- 
scale  atmospheric  tides  were  closely  coupled  to  the  dust 
content  of  the  atmosphere,  especially  during  great  storms 
(49:293).  The  Mariner  9  IRIS  (infrared  spectroscopy) 
experiment  showed  the  dust  is  about  60±10%  Si02  (18:430). 
Mariner  9  pictures  also  showed  dust  thrown  into  the 
atmosphere  as  high  as  40-60  km  (49:302).  Dust  storms  tend 
to  reduce  vertical  temperature  gradients  but  increase 
horizontal  temperature  gradients  (17:334),  to  warm  the  upper 
atmosphere  since  solar  radiation  is  absorbed  by  dust 
suspended  in  the  atmosphere,  to  increase  the  depth  and 
intensity  of  circulation  (17:335),  to  affect  the  polar  ice 
deposit  albedo,  and  to  vary  temperature  dependence  on 
latitude  (32:411).  The  following  figure  taken  from  (18:434) 
shows  some  of  these  effects: 
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Temperature  profile*  obi  niucd  by 
inversion  of  mtiimtmtnu  in  lit*  Cdi-cm”1 
COj  *b*orpiion  bond.  Tin  pronto*  were  obtained 
ftl  approximately  -30*  latitude  and  NOOhr 
Martian  local  time  during  tin  rcvolutioua  indi¬ 
cated.  During  revolution  30  tin  atmoapltera 
waa  heavily  dual  laden.  By  revolution  93  tin 
dual  waa  diaaipat  ing ;  clearing  waa  well  advanced 
by  revolution  174. 


Figure  1.1:  The  Effect  of  Dust  Storms  on  Atmospheric 
Temperature  and  Pressure 

During  dust  storms,  temperature  distribution  is  roughly 
constant  (29:389).  Furthermore,  during  the  decay  phase,  the 
atmosphere  is  more  isothermal  and  thus  more  stable  (49:302). 

The  lack  of  oceans  also  affects  the  atmosphere  in  ways 
not  originally  foreseen.  The  absence  of  seasonal  heat 
storage  by  oceans  makes  seasonal  temperature  change  very 
large.  Furthermore,  outgoing  infrared  radiation  is  not 
constant  with  latitude  but  decreases  from  equator  to  poles. 
Finally,  atmospheric  wave  character  is  changed,  resulting  in 
varying  circulation  regimes  (21:130). 

The  atmosphere  is  very  thin.  Average  surface  pressure 
is  less  than  0.01  that  of  earth.  However,  because  of  large 
variations  in  surface  topography  altitude,  some  locations 
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have  almost  rive  times  the  surface  pressure  of  others.  This 
fact  plus  meteorological  processes  make  the  surface  pressure 
highly  variable.  The  total  mass  of  the  atmosphere ,  unlike 
Earth's,  undergoes  major  seasonal  variations  of  as  much  as 
25-30%  as  C0<>  "snows  out"  at  the  winter  pole  and  is  released 
from  the  receding  summer  polar  cap.  Yet,  in  spite  of  its 
thinness,  the  atmosphere  still  supports  three  different 
kinds  of  clouds  as  well  as  dust  storms  ( 22:4280). 
Furthermore,  the  presence  of  dust  plays  an  important  role  in 
establishing  atmosphere  structure. 

The  Martian  atmosphere  is  usually  divided  into  two 
layers,  a  lower  neutral  layer  and  an  ionosphere  ( 21:228}* 
The  lower  atmosphere  consists  mainly  of  carbon  dioxide.  The 
characteristic  temperature  is  about  200  K  (21:28),  which  is 
generally  colder  than  the  average  surface  temperature  which 
is  about  250  K  in  the  daytime  (21:25).  This  difference  is 
due  to  the  presence  of  temperature  gradients  which  have  an 
average  value  of  about  1.6  K  per  vertical  kilometer 
(39:4363).  Higher  up  in  the  atmosphere,  large  daily  waves 
of  temperature  and  density  have  been  found,  indicating  a 
confused  and  variable  structure  (39:4364).  Below  125  km  in 
altitude,  where  the  atmosphere  is  well  mixed  by  turbulence, 
Viking  mass  spectrometers  have  discovered  that  96.5%  of  the 
atmosphere  by  weight  is  CO2  (34:4341).  The  composition  of 
the  lower  Martian  atmosphere  is  as  follows  (37:4637): 


9 


Tabic  1.2:  Composition 

of  the  Martian  Atmosphere 

( fraction 

by  weight) 

Carbon  Dioxide  { CO2  * : 

96. 5X 

Molecular  Nitrogen  (N2I: 

1.8X 

Arson  (Ar): 

1.5X 

Molecular  Oxygen  (02): 

• 

O 

Carbon  Monoxide  (CO):* 

0.05X 

Water  Vapor  UteOJ:0 

0.G2X 

Neon  (Ne): 

0.0001X 

Krypton  IKr): 

0.00003X 

Xenon  ( Xe  ) : 

0.00002X 

*  Uncertain,  variable,  or 

unevenly  distributed. 

Trace  amounts  of  other  gases  also  appear  due  to 
photochemical  reactions.  The  lower  atmosphere  supports  the 
ionosphere . 

The  topmost  layer,  or  ionosphere,  is  that  part  of  the 
atmosphere  where  constituent  densities  are  low  enough  for 
the  components  to  separate  diffusively  according  to  mass. 
Unlike  the  Earth’s  ionosphere,  its  peak  density  of  about  10° 
lons/cm  occurs  near  130  km  and  is  formed  primarily  of  O2 
ions  (approximately  90X  0g+  and  10X  C02+ )  that  have  been 
produced  in  local  chemical  reactions  (19:4351).  Another 
important  constituent  of  the  outer  limits  of  the  atmosphere 
is  atomic  hydrogen,  H,  which  is  detectable  over  10,000  km 
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above  rhe  surface  (37:210),  The  following  figures  taken 
from  l  I  *- 1 fe  1 J  provide  data  on  the  ionosphere: 


tUCKO*  MUMM4  StNUTT  u*’*| 


Water  Is  a  minor  constituent  of  the  atmosphere  due  to 
low  atmospheric  pressure  and  low  surface  temperatures, 
Nevertheless,  it  seems  to  play  an  important  role  in 
atmospheric  chemistry  and  meteorology.  The  atmosphere  is 
effectively  saturated  with  water  vapor,  yet  water  exists 
only  in  vapor  or  solid  form.  Viking  also  found  water  vapor 
is  uniformly  mixed  up  to  altitudes  of  10-15  km  and  shows 
strong  latitudinal  gradients  that  depend  on  season  and 
location  (13:4225).  A  strong  diurnal  variation  was  also 
detected  (13:4242).  The  total  amount  of  free  water  in  the 
atmosphere  is  estimated  to  be  equivalent  to  1.3  km  of  ice 
( 13:4246). 

The  vertical  structure  of  the  Martian  atmosphere,  i.e. 
the  relation  of  pressure  and  temperature  to  altitude,  is 
determined  by  a  complicated  interplay  among  several 
processes:  radiative,  convective  (natural  and  forced),  and 
advective  energy  transport  and  the  way  energy  from  the  sun 
is  at  first  absorbed  into  the  atmosphere  and  then  lost  by 
radiation  into  space  (5:161).  Because  of  the  wide  range  of 
density  in  the  atmosphere,  different  processes  are  dominant 
at  different  height  ranges.  Before  1965,  essentially  all 
knowledge  of  atmospheric  structure  was  based  on 
extrapolations  of  indirect  data  with  the  help  of  theoretical 
models  (21:25).  As  mentioned  previously,  many  of  these 
models  made  assumptions  that  proved  to  be  incorrect  and  so 
invalidated  the  extrapolated  data.  Space  probe  data  has 
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given  us  the  first  real  knowledge  of  the  structure  of  the 
Martian  atmosphere. 

Two  factors  control  the  vertical  structure  of  the  lower 
atmosphere:  its  composition  of  almost  pure  carbon  dioxide 
and  its  content  of  large  quantities  of  suspended  dust. 
Carbon  dioxide,  which  has  a  strong  infrared  spectrum, 
radiates  energy  with  great  efficiency  at  Martian 
temperatures  and  causes  the  atmosphere  to  respond  rapidly  to 
changes  in  the  amount  of  solar  radiation  received.  The 
suspended  dust  is  a  strong  absorber  of  heat  directly  from 
the  Sun's  radiation  and  provides  a  distibuted  energy  source 
throughout  the  lower  atmosphere  (29:386). 

Surface  temperatures  depend  on  latitude,  season,  and 
time  of  day  (39:4364).  (See  Fig.  1.4  Soviet  Mars  3 
photometer  data  (32:411)).  Viking  1  observed  a  wide 
temperature  fluctuation  from  day  to  night.  The  usual 
temperature  variation  was  from  about  150  K  just  before 
sunrise  to  a  high  of  240  K  in  the  early  afternoon  (42:3969). 
This  diurnal  variation  is  greatest  near  the  ground  and 
occurs  because  of  the  surface's  ability  to  radiate  its  heat 
to  space  quickly  during  the  night.  During  dust  storms,  this 
ability  is  impaired  and  the  daily  temperature  swing  is 
reduced  (26:484).  The  effect  of  dust  storms  on  pressure  and 
temperature  is  shown  in  the  Mariner  9  infrared  spectroscopy 
data  of  Fig.  1.5  (18:434). 
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Figure  1.4:  Atmospheric  Temperature  as  a  Function 
of  Time  and  Latitude 


Mariner  9  infrared  spectroscopy  data  also  showed  large 
diurnal  variations  { 15-30  K)  in  atmospheric  temperature  up 
to  about  30  km  and  winds  with  a  strong  tidal  component 
(18:423).  However,  other  oscillations  appear  throughout 
the  atmosphere  as  a  result  of  the  direct  input  of  energy 
from  the  Sun.  These  oscillations,  which  are  wavelike  in 
nature  and  are  akin  to  tides,  have  been  measured  as  pressure 
and  temperature  variations  and  lead  to  a  very  complex  wave- 
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like  rather  than  smooth,  veitical  structure  (  3 9 : -1 3 7 5 ) .  The 
mean  temperature  in  the  lower  atmosphere  falls  off  with 
altitude  at  a  lapse  rate  of  about  1.4  K  per  vertical 
kilometer  to  a  height  of  about  40  km  and  is  roughly  constant 
at  140  K  above  that  level.  This  low  lapse  rate  was 


unexpected  and  is  due  to  the 


Dependence  on  latitude  and  Martian 
local  time  or  atmoaplicric  temperature*  at  th* 
2-inlmr  level  (~  10  kin)  and  at  tltc  0,.1-mbar  level 
(-mikin).  The  temperature*  represent  average* 
over  data  obtained  on  revolution*  1-83  during 
the  diwt  Mtorm. 


large  amount  of  suspended  dust. 


LOCAL  TIMC  (MM.) 


Somo  a*  Fig.  9.  except  that  data  u*cd 
were  obtained  aa  tho  atmospheric  diut  \va* 
diaripating  on  revolution*  83-120. 


Figure  1.5:  Atmospheric  Temperature  as  a  Function 
of  Time,  Latitude,  and  Dust  Storm 

Mariner  9  radio  occultation  data  provided  sample  temperature 
gradients  (26:504): 
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Mxntpe  tcinj>cn\i«rc  gradient*.  o<l«l  revolution*. 

Figure  1.6:  Average  Atmospheric  Temperature  Gradients 

Above  about  100  km,  the  structure  of  the  atmosphere  is 
determined  by  the  onset  of  diffusive  separation  of  the 
various  atmospheric  molecules.  The  tendency  of  heavier 
molecules  to  lie  below  the  lighter  overcomes  the  tendency  of 
turbulence  to  mix  the  atmospheric  gases.  Absorption  of 
ultraviolet  light  from  the  sun  then  dissociates  and  ionizes 
the  constituents  and  leads  to  complex  sequences  of  chemical 
reactions  (5:27).  The  net  result  is  the  embedding  of  an 
ionized  layer  near*  130  km  and  an  increase  in  temperature, 
particularly  of  the  ionic  component,  at  higher  altitudes 
(19:4351).  The  top  of  the  atmosphere,  or  exosphere,  is 
characterized  by  a  temperature  of  about  325  K  averaged  over 
the  Martian  year  (44:469).  A  good  summary  of  the  ion 
reactions  is  provided  in  the  following  tables  (5:255  and 
5:275) : 
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Table  1.3:  Some  Ionosphere  Reactions 

SOME  IMPORTANT  CHARGE 
REARRANGEMENT  REACTIONS 


mo- 

4-  CO; 

-  CO 

4-0;* 

-f  1*22  eV 

(2)0* 

+  0; 

-0 

4-  o: 

4-  1*34  eV 

(3)0* 

+  N; 

-  NO* 

4-  N 

4*  N2cV 

(4)  Nr 

4-  CO. 

-  N. 

4-  CO, 

4-  1*79  eV 

(5)  Nf 

4*  0; 

-  N; 

4-  0. 

4-  3*50  eV 

(6)  N; 

+  0 

-  NO* 

4-  N‘ 

4-  3  08  eV 

(7)  CO* 

4-  CO, 

-  CO 

4-  CO  ; 

4-  0-22  eV 

(S)  CO- 

+  o. 

-  CO 

4-  0  ,* 

4-  1*93  cV 

19)  CO¬ 

4-  0; 

-CO, 

4-  or 

4*  1-71  cV 

HO)  A* 

+  CO; 

— *  A 

4-  CO  r 

4-  1*97  cV 

(II)  A’ 

+  CO 

-  A 

4-  CO* 

4-  1*75  eV 

(12)  A* 

4>  N; 

-  A 

+  N; 

4-  0*18  eV 

(13)  A* 

4*  0; 

-  A 

4-  0; 

4.  3  68  eV 

(14)  or 

4-  e 

-  0 

4-  0 

4-  6-96  cV 

(13)  CO* 

4-  e 

-  C 

4-  0 

4*  2-90  eV 

(16)  Nj 

4- 1 

-  N 

4-  N 

4*  5-82  cV 

(17)  NO* 

'  4 ■  < 

-  N 

4-  0 

4-  2*74  eV 

(18)  CO: 

1  4-  t 

-  CO 

4-  0 

4-  8*33  eV 

Table  1.4:  Some  Ionosphere  Constituents 


SOME  ION  STATES  AND  ASSOCIATED  THRESH- 
OLDS  FOR  PHOTOIONIZATION 


Parent 

Neutral 

Product 

Ion 

Electron 

Removed 

Threshold 

ineV 

0(JP) 

0*HS) 

2p 

13*6 

0*(:D) 

2p 

16-9 

0*('P) 

2p 

18*6 

0*(4P) 

2s 

28*5 

0+(JP) 

2s 

400 

A('S) 

A*!^) 

2p 

15*8 

A*(JS) 

3s 

29*2 

o  ,tJr;) 

0,*(-*/7.) 

•Vp 

12*1 

0;(V7.) 

•7.2  p 

1  A*  1 

0,(-77.) 

.7.2  p 

170 

o:(\r,-) 

<7.2  p 

18*2 

OiOffj 

7.2  s 

-28 

0.*(V/,) 

0.2  s 

-40 

NS('J;) 

N;(-r;> 

7.2  p 

15*6 

N«J//.) 

.7,2  p 

16*7 

N;(Jr;) 

0.2  s 

18  3 

7.2  s 

-35 

CO  ('IT*) 

co*<\r*) 

o2p 

140 

C0*(J//() 

.-Tip 

16-6 

C0*(2-*) 

o2s 

19-7 

CO;  (T;) 

co ;  c-//*) 

Irr. 

13-8 

C0:(-/7y) 

I.T. 

17*3 

co:(-r;) 

3o. 

18-1 
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All  atmospheric  motions  are  ultimately  driven  by 
pressure  forces  induced  by  thermal  gradients  and  supported 
by  complex  feedback  mechanisms  (5:110).  Hence,  atmospheric 
circulation  is  driven  by  heating  due  to  absorption  of  solar 
and  infrared  radiation  by  dust  and  CCb  as  well  as  sensible 
heat  exchange  with  the  ground  (17:322).  Furthermore, 
the  global  pattern  of  the  atmospheric  circulation  of  Mars 
shows  many  superficial  similarities  to  the  Earth,  although 
the  causes  of  motion  arise  from  different  circumstances. 

The  most  important  of  these  are  the  ability  of  the 
Martian  air  to  adjust  rapidly  to  local  conditions  of  solar 
heat  input,  the  lack  of  oceans  which  have  a  large  thermal 
inertia,  the  great  range  in  altitude  on  the  surface,  the 
strong  internal  heating  of  the  atmosphere  due  to  suspended 
dust,  and,  finally,  the  seasonal  deposition  and  release  of  a 
large  fraction  of  the  atmosphere  by  the  polar  caps  (49:298). 
Knowledge  of  winds  is  based  on  inferences  from  a  global 
series  of  measurements  made  by  Mariner  9.  These  data  can  be 
related  to  the  zonal  wind  field  if  we  assume  terrestrial 
conditions  (18:423).  The  predicted  wind  fields  show  a 
strong  dependence  on  the  Martian  season  since  large 
horizontal  temperature  gradients  build  up  across  the  edge  of 
the  newly  forming  polar  cap  in  autumn  and  winter  (49:298). 
The  result  is  the  formation  at  high  latitudes  of  strong  jet 
streams  with  eastward  velocities  in  excess  of  100  m/s 
(37:234).  On  Mars,  unlike  the  Earth,  there  is  also  a 
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relatively  strong  meridional  circulation  that  transports  the 
atmosphere  to  and  from  the  winter  and  summer  poles  I  5 : 1 2 i  I . 

The  general  circulation  pattern  is  occasionally 
•instable  and  exhibits  Large-scale  wave  motions  and 
instabilities  (5:1211-  A  wave  pattern  such  as  this  is 
cal  led  a  Rossby  regime  and  is  characterised  by  baroclmic 
instability,  which  has  been  observed  with  Mariner  9  TV 
observations  (29:373).  The  three  dimensional  wave  pattern 
is  shown  in  Fig.  1.7  (5:116).  Smaller  scale  wave  motions 
and  tides,  driven  by  the  Sun  and  surface  topography, 
include  a  strong  diurnal  component  and  are  present 
everywhere  in  the  atmosphere,  resulting  in  the  complex 
vertical  structure  (18:439).  This  turbulence  is  probably 
why  the  atmosphere  maintains  a  uniform  composition  to 
extremely  low  pressure  levels  and  is  presumably  also  an 
important  factor  in  maintaining  the  large  quantity  of  dust 
supported  in  the  atmosphere.  Occasionally,  for  reasons  not 
fully  understood,  the  atmospheric  dust  develops  into  planet¬ 
wide  global  dust  storms  that  obscure  the  planet  and  reduce 
atmospheric  transmission  (49:297). 

In  summary,  the  recent  space  probe  data  indicate  that 
the  atmosphere  of  Mars  is  both  position  dependent,  through 
latitude  and  longitude  effects,  and  time  dependent  through 
such  effects  as  these:  diurnal  and  seasonal,  annual  motion, 
solar  activity,  and  dust  storm  effects. 
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-Vertically  exploded  view  of  a  nave  of  short  length.  The  lower  portion 
shows  a  single  streamline  of  the  wave  while  the  upper  portion  shows  (1)  the  distri¬ 
bution  of  potential  temperature  in  two  vertical,  north-south  cross  sections  through 
the  inflection  points  of  the  wave  and  (2)  the  displacements  of  air  produced  by  the 
wave  fdouble  arrow  s).  In  a  short  wave  the  slope  of  the  displacements  is  greater  than 
the  slope  of  the  potential  isotherms. 


Figure  1.7:  Roaaby  Circulation  Regiae  in 
Three  Dimension* 
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ILL _ Ataosphere.-Models 


To  study  the  effects  of  aerodynamic  forces  on 
trajectories  at  orbital  speeds,  it  is  necessary  to  model  the 
planetary  atmosphere  in  which  the  flights  take  place.  The 
overriding  feature  of  the  atmosphere,  as  far  as  the 
spacecraft  is  concerned,  is  the  density  (-17:1).  Therefore, 
the  concern  in  modeling  the  atmosphere  is  to  conveniently 
and  accurately  represent  the  density. 

This  study  will  not  model  all  the  atmospheric  effects 
present  in  the  .Martian  atmosphere  which  contribute  to  its 
density.  It  is  meant  to  be  an  initial  attempt  at  modeling 
the  Martian  atmosphere  dynamically.  Only  long  term  effects 
are  kept;  short  term  fluctuations  are  averaged  out  or  an 
average  value  is  assigned  for  the  effect.  There  is  some 
doubt  in  much  of  the  work  on  Mars  due  to  assumptions  that 
are  not  proven  to  be  completely  accurate  as  yet  (5:18>l). 
As  a  result,  empirical  data  will  be  used  to  provide  typical 
or  reasonable  values  whenever  possible.  Both  dynamic  and 
static  models  use  1000  km  for  the  limit  of  atmospheric  drag 
effects;  while  the  ionosphere  has  hydrogen  ions  as  high 
above  the  surface  as  10,000  km  the  first  ions  that  cause 
significant  loss  of  momentum  upon  impact  are  located  at 
about  1000  km  (37:210).  Finally,  both  models  contain  some 
error  in  their  calculation  of  pressure  due  to  the  fact  that 
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pressure  is  supposed  to  be  calculated  with  respect  to 
altitude  from  the  local  Martian  surface,  not  with  respect  to 
distance  from  the  mass  center.  However,  more  detailed 
knowledge  of  topography  would  be  necessary  and  the 
complexity  is  beyond  the  intended  scope  nf  this  thesis, 
instead,  the  dynamic  model  uses  a  reference  pressure 
spheroid,  and  the  static  model  essentially  uses  the  Martian 
center  of  mass. 

Dynamic  Madel 

Because  the  mass  density  of  the  Martian  atmosphere  is 
heavily  time  dependent,  an  accurate  calculation  of  the 
atmospheric  forces  on  low  altitude  Martian  satellites 
requires  a  time  dependent  atmospheric  model.  This  highly 
variable  nature  of  the  Martian  atmosphere  means  predictions 
of  mass  density  at  a  specific  time  at  a  specific  position 
are  likely  to  be  greatly  in  error  (9:1).  However,  over  a 
long  term  period,  mass  density  can  be  predicted  with 
acceptable  accuracy. 

This  atmosphere  model  includes  seasonal  and  diurnal 
variation,  dust  storm  effects,  solar  activity  effects,  and 
the  effect  of  the  orbital  position  of  Mars.  The  change  in 
surface  pressure  due  to  polar  cap  sublimation/condensation 
is  viewed  as  averaged  out  over  the  long  term  of  these 
orbits.  Furthermore,  near  surface  variations  and  variations 
due  to  such  effects  as  topography  are  viewed  as  short  term 
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effects  only.  Finally,  the  complex  wave  perturbation  motion 
in  the  atmosphere  anti  related  latitude  and  iongituue 
gradient  representations  are  Left  for  follow-on  work. 

Reference  surlaee 

U1  altitudes  are  measured  with  respect  to  a  ••  l 
millibar  pressure  spheroid  (15*4321).  Since  Mars  has  no 
oceans  it  has  no  sea  lei  el  to  use  as  a  reference.  Instead, 
a  reference  pressure  spheroid  Is  determined  by  picking  the 
spheroid  which  In  a  Least  squares  sense  best  fits  Mariner 
occuLtation  data  for  a  6.1  mbar  isobar  surface.  This 
spheroid  nearly  coincides  with  the  average  surface,  ignoring 
topography,  and  is  specified  by  these  ellipsoidal  serainxes: 

A  =  3394.67  km 

B  =  3393.21  km  (3.1) 

C  =  3376.78  km  , 

where  4,B  are  in  the  equatorial  plane,  C  is  polar,  and  the  A 
axis  is  oriented  to  the  longitudinal  angle  108.5'  (19:324 

and  6:4353).  The  areocentric  coordinate  system  is  the  Mars- 
centered  coordinate  system;  it  is  equivalent  to  the  Earth- 
based  longitude  and  latitude  system.  The  fundamental  plane 
is  the  Martian  equatorial  plane,  and  the  principal  axis  is 
taken  as  pointing  towards  the  Martian  vernal  equinox 
(14:148).  The  Prime  Meridian  for  Mars  is  set  at  a  small 
crater  named  Airey-0.  The  celestial  longitude  of  the  Sun, 
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Ls.  in  areocentric  coordinates  is  an  angular  measure  of  the 
glance’s  orbital  position- 

Act  ivitv.  Variation  Effects 


The  Martian  atmosphere  is  extremely  sensitive  to  solar 
radiation  and  responds  quickly  to  any  changes  in  radiation 
intensity  (32:420).  Both  temperature  and  density  are 
affected  (36:214).  The  two  largest  effects,  the  l  i .  0*1  year 
sunspot  cycle  and  the  26.35  day  period  of  the  Sun’s  rotation 
(the  synodic  period  with  respect  to  Mars),  can  be  modeled 
for  an  average  intensity  over  a  long  term  period  19:3  and 
17:2). 

•V  common  measure  of  solar  radjation  intensity  used  in 
the  Literature  is  the  solar  F^q,?  centimeter  radio  flux, 
usually  given  In  units  of  10**^W7m^Hz  (19:4362).  Following 
the  method  in  reference  (10),  can  be  modeled  using 
empirical  data  to  create  an  average  amplitude  function: 
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where  PI=3 . 14 1592654 ,  Y  is  the  time  interval  in  years  from 
April  1,  1976,  d  is  the  time  interval  in  days  measured  from 
a  reference  date  chosen  to  match  available  Viking  1,2  Lander 
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data,  and  B  is  0.0011  times  the  value  of  with  the 

short-term  variation  averaged  out,  i.e.,  B  is  the  first 
bracketed  term  times  0.00U.  This  equation  is  suitable  for 
long  term  orbits  and  will  be  used  in  various  other 
equations . 

Effect  on  Temperature 

The  F^Qt7  solar  radiation  intensity  is  used  to 
calculate  the  exosphere  temperature,  T„ .  The  exosphere 
begins  where  the  moan  free  path  of  atmosphere  constituents 
equals  their  scale  height  (21:28).  This  calculation  is 
modeled  using  empirical  data  for  amplitudes  and  modeling  the 
distance  of  the  planet  from  the  sun  (8:3).  The  equation 
representation  is 


where  (a/R)  is  the  ratio  of  the  mean  solar  distance  to  the 
current  solar  distance,  and  the  radiation  intensity’s  decay 
with  the  square  of  the  distance  is  included  through  (a/Rp. 

Effect  on  Density 

Compared  to  the  Earth,  the  'an  atmosphere  is 
transparent  to  energy  radiation  and  witn. ut  oceans  does  not 
store  heat  as  readily.  As  a  result,  the  local  effect 
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differs  from  that  on  Earth  (-I2:39G9).  The  response  of  the 
Martian  atmosphere  to  solar  radiation  is  a  function  of  time 
and  can  be  broken  down  into  two  effects:  seasonal 
variations  and  diurnal  variations.  Furthermore,  time  can 
be  treated  as  a  function  of  the  local  angle  of  the  Sun. 

Seasonal  Variation 

A  simple  way  to  relate  this  effect  to  time  is  to 
express  it  as  a  function  of  the  difference  in  latitude  of 
the  Sun,  SNLAT,  and  the  latitude  of  the  spacecraft,  SCLAT. 
This  can  be  modeled  as  a  correction  to  the  turbopause 
altitude,  Zf,  where  Zf  is  defined  as  the  upper  limit  of  the 
thermosphere  or  lower  atmosphere  (8:3).  The  amplitude  of 
this  effect  is  very  small  {21:25  and  39:4367),  and  in  the 
absence  of  better  information  has  been  set  to  €s  =  0.1  km. 
The  empirical  equation  in  terms  of  turbopause  altitude  is 

dZf  =  €s  cos ( SNLAT  -  SCLAT).  (3.4) 
Diurnal  Variation 

The  diurnal  effect  is  expressed  as  a  function  of  local 
time  in  radians  by  relating  time  to  the  difference  in 
longitude  of  the  Sun,  SNLONG,  and  the  longitude  of  the 
spacecraft,  SCLONG.  In  addition,  (8:3)  adds  a  two  hour  lag, 
set  equal  to  PI=3.14...,  similar  to  that  seen  in  the  Earth’s 
atmosphere.  Thus,  the  resulting  equation  for  time  is 
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TIME  =  {SCLONG  -  SNLONG )  +  PI  . 


(3.5  1 


However,  this  effect  falls  off  with  latitude  until  it 
becomes  negligible  in  the  region  of  the  poles  since  less 
sunlight  reaches  there,  i.e.,  temperature  is  latitude 
dependent  (21:25).  As  a  result,  the  final  equation  will 
include  a  latitude  effect,  cos(SCLAT),  chosen  to  make  this 
effect  vanish  at  the  poles. 

Viking  lander  data  showed  the  diurnal  effect  decreasing 
with  decreasing  altitude,  being  only  a  minor  variation  at 
the  surface  (39:4367  and  47:2).  As  a  result,  this  effect  is 
deemed  negligible  in  the  lower  atmosphere  and  is  treated 
only  as  a  variation  in  the  upper  atmosphere.  There,  the 
effect  is  separated  into  two  components: 

1)  A  correction  to  turbopause  altitude,  Z^: 

As  with  the  seasonal  effect,  the  amplitude  of  this 
correction  is  small  (18:439)  and,  in  the  absence  of  better 
information,  is  set  to  €D  =  0.1  km.  According  to  (10:331), 
the  best  empirical  fit,  including  latitude  drop  off,  is 

AZf  =  €D  cos(TIME  -  3.6651914  )cos(SCLAT) .  (3.6) 

2)  A  variation  in  the  mixing  ratio  of  monatomic  oxygen 
above  the  thermosphere: 
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From  empirical  data,  monatomic  oxygen  is  the  dominant 
constituent  at  this  altitude  (34:4348).  This  variation  is 
accepted  in  (10:331  and  31:4384)  aa  being  represented  by  the 
following  equation  for  the  mixing  ratio  at  the  turbopause: 

f0  =  0.02  (1  -  C.5  sin(TIME)/cos(SCLAT) ) .  (3.7) 
Dust  Storm  Effects 

While  dust  is  always  present  in  the  Martian  atmosphere, 
the  primary  turbidity  effects  are  those  of  the  global  duat 
storms  which  occur  annually  as  Mars  passes  perihelion 
(38:262).  Global  dust  storms  are  those  dust  storms  that 
grow  out  of  local  dust  storms  at  certain  preferred  locations 
and  seasons  to  obscure  the  entire  planet  and  throw  dust  into 
the  atmosphere  as  high  as  60  km  (49:302).  They  cause  both 
density  and  temperature  changes.  However,  observed 
variation  from  year  to  year  has  been  great  and  can  be 
expected  to  remain  so  (49:296).  Therefore,  the  only  way  to 
model  this  effect  is  to  choose  a  typical  year  global  dust 
storm  as  the  "average"  effect.  For  this  purpose,  Viking 
data  on  the  1977  storms  is  chosen.  These  storms  showed  a 
double  peak  with  maximum  thicknesses  of  7.5  km  with  onset  at 
solar  longitude  of  Ls  =  210'  and  13.5  km  at  solar  longitude 
of  Ls  =  270*  in  the  areocentric  system  (49:290).  Note,  the 
perihelion  of  Mars  corresponds  to  an  areocentric  solar 
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Longitude  of  about  Ls  =  250‘.  The  onset  of  all  global  dust 
storms  was  within  60'  of  this  longitude  (49:299). 

The  dust  storm  effect  appears  as  an  isothermal  layer 
inserted  at  the  bottom  of  the  Martian  atmosphere,  with  the 
higher  thermal  layers  pushed  upward  by  this  amount  (26:493). 
Therefore,  using  Viking  data  the  turbidity  layer  is  modeled 
as  (10:332) 


A  Z0S  =  1 . 2  kra+  2 


i  =  l,2 


AZi 

' 

0 . 65 

where  A  - 
AZ2  = 
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(3.9a) 
!  3 . 9b ) 
(3.9c) 
( 3 . 9d  } 


The  equation  uses  1.2  km  as  an  average  impact  on  the 
turbopause,  and  the  right  term  fits  the  empirical  data  for 
the  additional  overlapping  effects  of  the  two  peaks  at  210* 
and  270*  (4S:290). 


Annual  Motion  Effects 

The  annual  motion  of  Mars  about  the  Sun  changes  the 
distance  between  the  two  bodies  as  well  as  the  orientation 
of  the  planet  to  the  Sun’s  radiation.  These  changes  result 
in  variation  in  the  temperature  profile  of  the  atmosphere  in 
addition  to  changes  in  the  "sire"  (and  so  pressure)  of  the 
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atmosphere  —  the  atmosphere  alternately  swells  and  shrinks. 
However,  the  pressure  variation  is  less  important  than  the 
other  effects.  Annual  motion  is  basically  a  time  effect 
which  can  be  expressed  as  a  function  of  both  the  longitude 
of  the  Sun  in  areocentric  coordinates,  Ls,  (measured  from 
the  Martian  autumnal  equinox,  which  corresponds  to  a 
heliocentric  longitude  for  Mars  of  85*),  and  the  Mars-Sun 
distance  R  expressed  nondimensionally  as  (a/R)  (8:4).  The 

dust  storm  effect  in  the  last  section  i3  basically  an  annual 
motion  effect. 

Empirical  data  show  a  good  approximation  for  the  annual 
motion  effect  on  pressure  is  (9:5) 

P  =  P0  (1  -  0.24  sin(Ls  -  59')).  (3.10) 

This  pressure  effect  can  be  included  by  correcting  the  6.1. 
mbar  ellipsoid  used  as  the  reference  level  in  this  model 
atmosphere  by  (10:332) 

AR6.1  s  -11.3  ln[  1  -  0.24  sin(Ls-59*)]  km.  (3.11a) 

In  recognition  of  the  uncertainty  in  this  term,  and  to  avoid 
repeated  evaluation  of  the  log  term,  this  is  approximated  to 
the  first  order  of  the  logarithmic  series  for 
ln(l-x)=  -(x+x^/2+. . . )=  -{ -11 . 3  )  ( -0. 24sin( ) )  by 

ARe.l  =  "2'8  sin(Ls  -  59')  km.  (3.11b) 
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This  is  based  on  an  average  temperature  over  the  entire  day 
at  the  reference  ellipsoid  of  220  K  (21:25). 

The  two  most  important  effects,  the  temperature  profile 
variation  and  the  change  in  size  of  the  atmosphere  are 
modeled  using  empirical  data.  Basically,  accepted  usage  is 
to  form  a  temperature  function  which  is  piecewise  linear 
between  specified  data  points.  One  such  function  is  given 
in  the  following  table  (39:4364  and  9:5): 


Table  3.1:  Piecewise  Linear  Temperature  Function 


of  the  Martian 

Thermosphere 

(Where  the  ratio  (a/R)  is  the 

semimajor  axis  of  Mars' 

orbit  divided  by  the  current 

Mars  to  Sun  distance. ) 

Altitude  (km) 

Temperature  (K) 

Reference,  R61MB,Z=0 

220* (a/R) 

ZL  =  A  ZDS( 2.75) 

220* (a/R) 

ZA  =  ZL.  +  32  (a/R) 

140* (a/R) 

ZB  =  AZds  +  100  (  a/R }  140  *  (  a/R) 


ZF  =  AZds  +  125(a/R)  +  A  ZSN  172.5*(a/R) 


where  A  Zpg  is  the  dust  storm  correction  and  A  Zg^  is  the 
seasonal  and  diurnal  correction.  The  table  represents  two 
isothermal  layers  and  two  layers  of  constant  temperature 
lapse  rates  (24).  Furthermore,  both  the  exosphere 
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temperature  and  the  solar  radiation  flux,  o . 7 »  depend  on 
the  annual  motion. 

Pena  Uji. .Calculations 

For  mass  density  calculations,  the  atmosphere  can  be 
divided  into  two  regimes:  a  lower  neutral  atmosphere,  called 
the  thermosphere,  and  an  upper  layer,  called  the  ionosphere 
(21:228).  The  dividing  line  between  these  two  layers  is 
called  the  turbopause  and  is  stated  as  an  altitude,  Zf. 
This  altitude  varies  between  approximately  115  km  to  150  km, 
with  an  average  value  of  about  130  km  (34:4346).  This 
average  altitude  also  represents  the  region  of  maximum  ion 
density,  i.e.,  this  is  where  a  spacecraft  would  "enter"  the 
atmosphere.  Consequently,  satellites  are  placed  in  orbit 
above  this  level  and  most  of  the  satellite  atmospheric  drag 
is  in  the  upper  atmosphere  regime.  The  two  regions  do  not 
behave  the  same  (18:435). 

The  lower  atmosphere  or  thermosphere  can  be  treated  as 
a  neutral  gas  with  an  ideal  gas  molecular  weight  of  u  =  43.3 
gm/mole  br.sed  on  Viking  ?,  data  (39:4368).  The  thermal 
layers  of  this  lower  atmosphere  vary  with  annual  motion  and 
a  temperature-altitude  model  is  given  in  Table  3.1.  Mass 
density  is  calculated  using  pressure  and  temperature 
information  in  the  ideal  gas  law.  The  pressure  in  this 
regime  is  known  to  follow  an  exponential  law  (5:232),  and 


32 


the  pressure  is  found  by  solving  the  barometric  lav 
analytically.  In  the  two  isothermal  layers,  the  pressure  is 
calculated  using  the  equation  (21:27) 


( ( -2-Zi }/H i ) 

P  =  PL  e, 


(3.12) 


where  P^  is  the  pressure  at  the  lower  bound  of  the  layer, 
Zjl  »  and  the  scale  height  for  the  layer  is  given  by 


Hi  =  kTi/ogi,  (3.13) 

where  k  =  Universal  gas  constant, 

Ti=  temperature  in  the  layer  (K), 
aii=  mean  molecular  weight 
gi=  local  gravity. 


In  the  two  layers  of  constant  temperature  lapse  rate  (linear 
temperature  function),  the  pressure  is  calculated 
differently.  The  generalized  equation  for  the  molecular 
weight,  m,  of  nonisothermal  regimes  changes  equation  (3.13) 
to  (39:4368) 


A  =  - 

k  T 

d£ 

P  g 

dz 

where  P  is  the  pressure  and  z  is  the  altitude.  Rearranging 
this  equation  and  multiplying  by  (dT/dT)  yields: 
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d£  =  (ag/kl  dz  dl  =  (Kg/k)  dz  dl 

P  T  dT  dT  T 


=  (ng/kT* )  dl,  (3.15a{ 

T 

where  T'  is  (dT/dz).  Integration  yields: 

lntP^/Pi)  =  (mg/kT’)  ln(T2/T1) 

(ag/kT' ) 

=  ln(T2/T1).  (3.15b) 

Taking  the  antilogarithm  produces: 

(agi/kT* ) 

P  =  PilT/Ti),  (3.16) 


where  i  refers  to  the  lower  bound  of  the  layer  as  before, 
and  T*  is  approximated  as  the  temperature  lapse  rate 
(AT/4  Z).  Once  the  pressure  is  known,  the  mass  density  can 
be  calculated.  The  standard  relation  between  pressure  and 
the  mass  density,  p,  is  (34:4344): 

dp/dz  =  -  pg .  (3.17) 

Substituting  this  into  equation  (3.14)  yields  the  equation: 

nP 

;i  _  — .  (3.18) 

kT 
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The  upper  atmosphere  or  ionosphere  does  not  follow  a 
simple  exponential,  ideal  law  formulation.  This  regime  is 
composed  of  dissociated  ions  and  molecules  which  are  not 
well  mixed  (19:4351).  Therefore,  the  mass  density  must  be 
calculated  using  properties  of  the  constituents  of  this 
layer.  The  most  straightforward  way  to  calculate  the  mass 
density  is  to  use  the  number  densities  of  the  constituents, 


modifying  the  process  with  a  measure  of  the  mixture  ratio 
between  them.  There  are  nine  constituents  included  in  this 
model  (see  Table  3.2). 


Tabic  3.2:  Atmosphere 

Model  Ionosphere  Constituents 

Constituent 

Molecular  Weight  (gm/mole) 

Carbon  Dioxide  (CO2) 

44.01 

Monatomic  Oxygen  (0) 

16.0 

Molecular  Nitrogen  ( N2 ) 

28.012 

Argon  (Ar) 

39.948 

Carbon  Monoxide  (CO) 

28.01 

Molecular  Oxygen  (O2) 

32.0 

Helium  (He) 

4.0 

Monatomic  Hydrogen  (H) 

1.0 

Molecular  Hydrogen  ( H2 ) 

2.0 

The  mass  density  is  calculated  in  this  region  using  the 
methods  found  in  references  (9)  and  (10).  The  number 
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densities  are  calculated  as  a  function  of  the  geopotential 
altitude,  x.  The  geopotential  altitude  is  defined  as: 

( 2 -zf J 

(3.U) 

where  z  =  altitude  of  spicecraft  above  pressure  reference 
surface 

2 j.*-  altitude  of  turbopause 

Rf =  areocentric  distance  of  turbopause. 

The  temperature  is  calculated  using 

-  X 

a 

T  =  T.  -  (T.  -  Tf)  e  (3.20) 

where  a  =  (T.  -  Tf)  /  (AT/AZ)f  =  T0/18,  (3.21) 

Tf  =  temperature  at  turbopause, 

(AT/AZ)f  “  approximation  for  T*  at  turbopause. 

The  partial  pressure  for  a  constituent  i  is  calculated  using 

“  X  A 

Hi  Hi 

pi  =  pf  fi  e  (Tf/T),  (3.22) 

where  fi  is  the  mixing  ratio  of  the  constituent  at  the 
turbopause,  Zf,  (34:4345)  and  Hi  is  the  constituent  scale 
height  defined  by  (34:4346) 
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Hi  s  kJL  , 


(3.23) 


with  an  being  the  molecular  weight  of  the  constituent  and  gf 
being  referenced  to  the  turbopause.  Now,  use  the  partial 
pressure  to  calculate  the  number  density  in  the  rewritten 
gas  law  equation: 

Pi  N 

ni  = -  (3.24  ) 

k  T 

Finally,  the  mass  density  is  found  by  summing  the  number 
densities  and  converting  to  mass  units: 

-  ni  “i 

U  =  i - .  (3.25) 

N 

where  N  =  6.023  x  i026  /kg  mole  ,  k  =  8314.39  J/kg  mol  K, 
and  i  =  for  the  nine  constituents.  The  mass 

density,  u,  can  now  be  used  in  drag  equations  to  determine 
atmospheric  drag  on  a  spacecraft  in  the  orbital  equations  of 
motion. 

Static  Model 

The  static  atmosphere  model  is  a  much  simpler  model  of 
the  atmosphere.  It  is  not  meant  to  be  a  highly  accurate 
representation  of  the  actual  Martian  atmosphere  but  rather  a 
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quick  means  to  obtain  an  upper  bound  approsimat ion  ot  the 
mass  ‘tensity.  A  number  ot  fundamental  assumptions  are  made 
following  the  methods  Antroctyceri  in  reference  l  >17).  These 
assumptions  great Ly  streamline  the  necessary  calculations, 
which  saves  valuable  computer  time  if  great  precision  to 
real  Martian  atmospheric  conditions  is  not  required. 

The  static  atmosphere  is  modeled  quite  differently  from 
the  dynamic  model.  The  effect  of  composition  of  the 
atmosphere  is  assumed  to  be  negligible,  although  that 
composition  is  very  important  in  understanding  the 
ionosphere.  Instead,  the  atmosphere  is  treated  as  a  uniform 
gas  of  unvarying  composition.  This  is  a  reasonable 
assumption  if  the  atmosphere  is  modeled  to  provide  an  upper 
bound  on  density  rather  than  match  varying  actual 
conditions.  Also,  a  simple  "spherical  symmetry"  modified  by 
elLipticity  is  assumed,  L.e.,  the  atmospheric  pressure  is 
essentially  a  function  only  of  radial  distance  from  the 
center  of  the  planet.  As  mentioned  before,  this  should  be 
a  function  of  altitude,  which  is  a  function  of  topography. 
However,  this  model  is  meant  to  be  a  very  convenient  way  to 
estimate  density.  Also,  solar  radiation  effects  on  the 
atmosphere  are  regarded  as  being  negligible.  Again,  this  is 
acceptable  within  the  framework  of  the  intent  of  the  model. 
Furthermore,  the  atmosphere  is  assumed  to  be  nonrotating. 
This  is  a  reasonable  assumption  since  the  aerodynamic  forces 
generated  by  the  rotation  have  larger  effects  only  at  low 
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altitudes.  One  measure  of  the  effect  is  to  examine  maximum 
rotational  speed  of  Mars  at  the  equator.  This  can  be 
represented  as  the  ratio  of  the  surface  rotational  speed 
divided  by  the  circular  orbital  speed  at  the  surface 
expressed  as  a  percentage.  For  Mars,  this  measure  is  given 
as  seven  percent,  which  is  a  maximum  (47:3).  Therefore, 
this  is  a  reasonable  assumption.  The  net  effect  of  the 
above  assumptions  is  that  there  is  r  *ime  dependence  built 
into  the  model.  Another  very  important  assumption  is  that 
the  density  decreases  exponentially  with  altitude. 
Combining  equations  (3.17  i  and  (3.18)  and  reorganizing  after 
taking  the  resulting  derivative  yields 


du  =  -  am  +  1  ill 

p  kT  T  dr 


dr 


(3.26a) 


=  -  0  dr, 


(3.26b) 


where  0  is  defined  as  the  bracketed  term  above  and  is  the 
reciprocal  of  the  scale  height  (1/H).  If  0  is  considered  to 
be  constant  over  some  small  altitude  interval,  the 
integrated  density  function  is 

-0( z-z0) 

p  =  Pq  e  (3.27) 

-( z-z0 )/H 

or  P  =  Po  e. 
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\s  a  result,  the  only  elements  of  the  model  that  must  be 
determined  under  these  assumptions  are  the  input  values  for 
the  reference  density  p0,  the  corresponding  reference 
aLtitude,  and  the  scale  height  H  for  the  altitude  interval. 
Unfortunately i  there  are  many  models  and  various  data  in  the 
literature  and  they  don’t  always  agree  closely.  Therefore, 
a  representative  set  must  be  chosen  from  empirical  data. 

The  scale  height  has  been  observed  to  increase  with 
altitude  from  about  10  km  at  the  surface  to  a  maximum 
average  ionospheric  scale  height  of  38.5  km  (26:484).  Since 
the  spacecraft  will  be  spending  its  atmospheric  time  in  the 
ionosphere,  the  scale  height  should  be  somewhere  near  this 
upper  number.  Reference  (5:156)  uses  Mariner  data  to  give 
an  upper  bound  value  for  scale  height  of  about  35  km. 
Therefore,  a  reasonable  scale  height  datum  in  this  range  is 
36  km  for  the  ionosphere.  According  to  Viking  1  empirical 
data  this  corresponds  to  a  reference  height  of  about  361  km 
and  reference  density  of  approximately  6  x  10-13  kg/m3 
(39:4373).  These  values  result  in  a  density  equation  of 

.  „  _  (361  kra  -  z)/36  km 
P  =  (6  x  10'iJ  kg/ra*1 )  e,  (3.28) 

This  is  the  equation  used  in  the  static  atmosphere 
model  to  calculate  the  mass  density. 
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JS-u  Orbital  Equation  Model 


A  model  for  orbital  motion  about  Mars  is  developed 
using  Keplerian  motion  (Lagrange's  Planetary  Equations) 
subject  to  certain  perturbations.  The  perturbations 
included  are  the  aspherical  gravity  field  of  Mars 
(geopotential  or  central  body  disturbing  potential),  the 
gravitational  attraction  of  the  sun  (third  body  disturbing 
potential),  atmospheric  drag  (using  two  different  models  to 
generate  atmospheric  mass  density  for  the  calculation),  and 
solar  radiation  pressure.  This  section  addresses  each  part 
of  the  overall  model  in  order  to  build  the  complete  thesis 
model . 


Lagrange's  Planetary  Equations 

Lagrange's  Planetary  Equations  are  usually  written  in 
terms  of  the  classic  orbital  elements  a,e,i,Q,w,  and  M. 
However,  there  are  disadvantages  to  using  these  elements  as 
the  variables  in  the  orbital  differential  equations.  The 
typical  Lagrange's  Planetary  Equations  given  in  (20:29)  are: 


da  _ 

2  ^  R* 

dt 

na  ^  M 

de  _ 

1-e2 

dR* 

_  /1-e2  ^  R 

dt 

na2e 

d  M 

na2e  h  e 

(4. i.la) 


(4.1.1b) 
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aw  _ 

cosi 

b  R 

+  v  l-e"  aR 

I  4  . 1 .  lc  ) 

dt 

na"sini*‘ l-e“ 

^  i 

na“e  h e 

d  r  = 

cos  i 

1 

^R* 

(1.1. Id) 

dt 

\*  *) 
na“  sinivl-e" 

^  w 

2  9 

na  siniv l-e“ 

Z  12 

d*J  = 

1 

$  R* 

( 4. L. lei 

dt 

2  ,  ,  2 
na  sini  , l-e 

Ti 

dM  = 

n  .  l-e2  ^  R*  - 

2 

c)  R* » 

( 4 . 1 .  If) 

dt 

na2e  ^  e 

na 

<)  a 

where  R*  is  the  disturbing  potential,  assumed  a  function  of 
the  classical  elements,  and  n  is  the  mean  motion. 

As  can  be  seen  above,  i  and  e  are  in  the  denominators 
of  some  of  these  differential  equations  and  will  produce 
singularities  if  they  become  vanishingly  small.  To  avoid 
this  problem,  two  actions  are  taken.  First,  to  avoid  the 
singularity  in  inclination,  the  inclination  is  never  set 
equal  to  zero.  Second,  new  orbital  elements,  designated  h 
and  k,  are  introduced  which  will  eliminate  the  singularity 
due  to  zero  eccentricity  (28:3-1): 

h  =  e  sinw "  (4.1.2a) 

k  =  e  cos w.  (4.1.2b) 

To  facilitate  the  numerical  integration,  a  singly 
averaged  method  is  used.  This  method  removes  the  fast 
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v.iptabLe  (terms  involving  the  mean  anomaly  MJ  in  the 
geopotential  and  third  body  potential  before  numerical 
integration.  However,  terms  involving  the  mean  anomaly  but 
giving  rise  to  resonant  terms  with  the  geopotential 
tesserals  are  retained.  This  is  exactly  what  we  want  since 
only  long  term  orbital  behavior  is  of  interest.  To 
emphasize  resonance,  the  mean  anomaly  will  be  replaced  with 
another  new  orbital  parameter  called  the  stroboscopic  mean 
node,  L$  (16:167-189). 

While  affecting  all  Keplerian  elements,  the  greatest 
effect  of  resonance  shows  up  in  the  longitude  dependent  L# 
(20:19-56).  The  stroboscopic  mean  node  is  defined  as 
follows : 


J±-±_ H 
s0 


(-1.1.3) 


where  Sq  is  the  ratio  of  two  relatively  prime  integers,  P/Q, 
approximating  the  number  of  nodal  crossings  per  planet 
revolution,  and  8g  is  the  Prime  Meridian  angle.  Reference 
(20)  provides  a  good  description  of  the  physical  meaning  of 
the  stroboscopic  mean  node. 

Numerical  quadrature  is  used  to  average  the  potential 
due  to  solar  radiation  pressure  and  to  drag.  The  resulting 
averaged  equations  of  motion  are  integrated  numerically 
using  a  multistep,  variable  order  and  variable  step-size 
integrator  (28).  Due  to  the  elimination  of  the  fast 
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variable  in  the  equations  of  motion,  large  step  sizes  of  the 
order  of  days  are  possible.  This  suits  the  long  term  orbit 
studies  desired,  but  eliminates  short  term  information  since 
that  information  has  been  averaged  out  (28:1-2}, 

The  orbital  equation  model  uses  a  planet  equator  of 
epoch  coordinate  system,  which  is  a  non-rotating  coordinate 
system  where  the  X-Y  plane  is  the  equator  plane,  defined  at 
some  epoch,  of  the  planet  being  orbited  by  the  spacecraft. 
Polar  motion  is  ignored.  Reference  {28}  provides  a  set  of 
recommended  orbital  elements  for  the  sun  with  respect  to 
Mars'  equator  and  equinox  of  epoch  and  the  location  of  the 
Prime  Meridian  for  three  epochs.  The  closest  epoch  to  the 
mission  is  chosen;  the  positional  error  introduced  is  less 
than  0.1‘  if  the  solar  ephemeris  reference  date  is  very 
close  to  the  mission  date  {28:2-1,2-11). 

Lagrange's  Planetary  Equations  must  be  transformed  from 
the  classical  Keplerian  orbital  element  set  to  the  new 
orbital  parameter  set  of  a,h,i,k,Q,  L jj.  Substituting 
derivatf  es  of  the  disturbing  potential  into  Equs.  (4.1.1) 
and  using  the  chain  rule  to  introduce  the  new  orbital 
elements  yields  (30:3-2): 

da  =  2  1  ^  R  (4.1.4a) 
d  t  na  Sq  d  Ltf 

dh  -  v"l-e2  ^R-  k  coti  dB_hB’/l-e2  >R  (4.1.4b) 
dt  na"  ^k  na2Vl-e&  <)  i  na2Sp  ^ 
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k  $R  -  h  d*  +  l  c>R 
^  h  ^  k  S  q  ^  £\j 

^  R  +  ^  R  (4.1.4c) 

^  Q  ^  £n 

dk  =  _  yl-e2  ^  R  +  h  coti  ^  R  _  k  P'/l-e2  ^  R  (4.1..  Id) 

dt  na“  na2y 1-e2  ^i  naSQ 

d£  -  1  ^  R  (  4  . 1 . 4  e  ) 

„  __  -  g  \ ~ 

dt  na~  sini  vi-e  0  i 


d£N  -  n  -  dOg  +  1  p 8 yi-e2  h  oR  r  k  «J  R 

dt  Sq  dt  na2S(j  ^  h  c>  k 


-  2a  ^  R  + 

Sq  -  cosi 

^  R  , 

(  4  . 1 . 4  f  ) 

^  a 

sini  /1-e2 

<)  i  ■ 

where  e  =  vh2  +  k2 , 

O 

n  =  /p/a  ,  and 

P’  = 

1 

1  +  /1-e2 


di  _  coti 

’  v  2 
dt  na“yl-e 


This  model  for  the  aspherical  gravity  field  perturbation 
(geopotential)  is  developed  in  detail  in  references 
(11), (20),  and  (48).  The  geopotential  is  expressed  by  the 
following  relationship  where  the  equation  is  made 
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dimensionless  using 
planetary  radius  Re: 


V  (  r ,  U ,  o )  = 


-n  ®  i  r 

_  V  V 

r  1=0  raSO 


the  gravitational  parameter  p  and 


r 

R, 


P^t  sinO  J  (C  jracosmo  +  S^smmo), 

U  .  2 . 1 J 


This  equation  is  transformed  into  a  Keplerian  variable 
form  using  two  new  functions  that  depend  on  Keplerian 
elements:  the  inclination  function  Fjmp(i)  and  the 
eccentricity  function  Gjfpqte}.  The  complete  derivation  may 
be  found  in  (11)  or  (20). 

The  inclination  function  is  defined  as  follows: 


„  tmax 

F  1  =  £ 

*  imp1 - •  t=0 


( 21-2t ) ! 


t!(l-t)!(l-m-2t)!22i“2t 


m 

JC 

s=0 


m 

s 


cos s(i)  I 
c 


l-m-2t+s 

m-s 

c 

p-t-c 

Sinl-m-2t( ^  j 


(-l)c"k,  (4.2.2) 


where  c  =  all  integers  where  binomial  coefficients  are 
nonzero.  The  restrictions  on  binomial  coefficients  result 
in  a  lower  limit  on  index  c  of  ci=max(  0 ,  p-t-m-s ) ,  and  an 
upper  limit  of  cu=min( 2-m-2t+s , p-t ) ,  cu  >  c^. 

The  eccentricity  function  is  defined  as  follows: 


Glp(2p-1)= 


1  p’-l 

i^+l ( l-e2 )1~l  d=° 


r 

i-i 

2d+l-2p ’ 

2d+l-2p ' 

d 

2d+l-2p’ 

(4.2.3) 


46 


where  p' 


p  for  p<  1/2 

i-p  for  p>  i/2 , 


Now  the  disturbing  function  may  be  expressed  as 


•  1 

R  =  I  V1m, 

1=2  m=0 

(4.2.4) 

whe  re 

v*1 

vinraIn  p!o^^mP^q=-«  [Glpq^e  *Glmpq(  wiM,Q,Q  )j 

(4.2.5) 

and  simpq= 

C|mcosmo  +  Sjmsinmof  i-m  even 
-Sjmcosmo  +  Cjmsinmo,  i-m  odd 

(4.2.6) 

with  o  =  ( 1- 

-2p)w  +  (l-2p+q)M  +  m(Q-0). 

(4.2.7) 

Alternately,  converting  to  the  modified  elements  using  the 
equations  and  definitions  in  (28:3-4,7),  V^m  becomes 

dlmp?0Flrap(  ^q=5«§  '  lpql  e  ^Imq^1*  hi  ^  j  *  (4.2.8) 


Now  substitute  the  derivatives  of  Vjlb  with  respect  to 
the  modified  orbital  element  set  into  Lagrange’s  Planetary 
Equations  (4.1.4)  to  yield: 


da 

dt 


2 

naSg 


1 

»Re  1 

71*1  "pio 

cl 


0  ,  m  =  0 . 


1 

Flmp  G*lpq  ” 
at 


h  ^imq 
^  lN  ' 


m  not  0 


(4.2.9a) 
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or  i 

vl-e“  pRe 


2  xni  |=oFi«p  r'^pq  ^  k  + 

•  l* 


}simq  k  dG' 

V  +  simq 

d  k  e  de 


k  coti  PRe  1  ^  Rlmp 

i  ii  i  1. 1  ■  ■  .  .  t  *  v  -  -  -  n  > 


2  * ,  2  1  +  1  J-*mpIo  \  . 

v  1-e  a*T*  v  u  d  l 


G  lpq^e^imq 


hP 1 y 1-e2 f  pRg  i 


»PM-^  pRe  1  . 

na2S0  a-*+1  iW  ^  FimP  °  lp<1  m 


0  for  m  =  0 


(4.2.9b) 


di  coti  k  pRe  J la  l  ^  <^Sjmq  ^  h  dC 

=  ~2~  9  i+T  “7  '  pio  FlmP  G’iP9  TT  *  S**<1  ~  T~ 

dt  na  y l-e  a  sini  H  L  d h  e  de 


coti  h  HR,  J i„  J  ,  zslmq  k  do1 

n*2A-c2  .*«  sini  »?°  l  Up  lpq  <>  k  +  Sj"'1  I  do 


coti  1 

S(jna2sini/l-e2  na2  sini  /1-e2 


pRe  i  1  <5simq»  m  not  0 

x  “I+T  Ji»  mpIo  Fi®p  G'ipq  “  77 

aA  *  v  m 

0 ,  m  =  0. 


(4.2.9c) 


o  1 
dk  -  /l-e^  pRe 


e  i  „  <*ai*q 

2  1+1  Jim  p=0  FlraP  °  XP9  \.  +  Si®9 


na  a 


h  coti  pRe  ^  j  aFlmp  _ 

+  2  *.  2  1+1  Jim  P=0  A  •  G  -zP<I(e,S/Bl<l 

na  y  l-e&  a  F  ox 


h  dG’ 


e  de 
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-  k  |3\L-e2  pRg 


nnSg 


,1+1  J-*m  mPlo 


d!3  _ 


pRg 


dt 

V 

na“ 

sin 2  v 1 

-e2  a 

d% 

_  n 

-  d0g  + 

l 

dt 

s0 

dt 

na2SQ 

r  r 

Flmp  G ' lpq  “ 
m 

0,  m  ■ 

=  0. 

J^®PIo 

£Flmp 

^2 

imq 1 


LN 


( -1 . 2 . 9d  } 


( 4 . 2 . 9e  J 


ON'l-e2 


x 


2 

URe 


,i+L 


P=0 


^G2mq  ^  dG’ 

x"p  lM  dh  *  Sj"“  ;  ^ 


+  k 


a 


2 

PRC  j 

2+1  J2«p?o 


_  ^S2»q  k  dG’ 

Flmp  G  lpq  v  +  s2mq 

dk  e  de 


r  2 

2 

+  2a 

2+1  l2+2Jimpio 

Flmp( A  ^G  * lpq ^  e )S2mq 

>0  “  cosi  pRe 


■  •  2  2  +  1  J^ranIo 

sini  vl-e  aA  A  p  u 


^F1 


mp 


^  i 


G'lpq(e)s2mq 


.  { 4 . 2 . 9f ) 
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Solnr  Radiation  Pressure  Model 


The  acceleration  of  a  satellite  due  to  direct  solar 
radiation  pressure  is  given  by  { 28 : 3-29 J : 


A’ 

r  =  r  ps  - 


m 


sun 


■  v» 


(4.3.1) 


where  £=  satellite  position  vector  relative  to  planet  center 
rsun  =  distance  of  the  sun  at  1  A.U. 

Es=  position  vector  of  the  sun  relative  to  planet  center 
JCv=  i :  -  £s 

ps=  radiation  pressure 

t  =  dimensionless  black  body  radiation  constant 
A'  =  effective  spacecraft  area  for  radiation  pressure 
m  =  mass  of  the  satellite 
r  =  unit  vector 


For  an  orbiting  spacecraft,  a  standard  approximation  is 
to  take  rv  *  -  £s.  Also,  since  the  orbital  eccentricity  of 
Mars  is  small,  0.093387,  assume  rs  =  ap,  where  ap  is  the 
semimajor  axis  of  Mars’  orbit.  Further  assume  rsun  =  ae» 
where  ae  is  the  semimajor  axis  of  the  Earth’s  orbit,  since 
its  eccentricity  is  0.0167.  The  effect  of  these 
assumptions  is  to  change  equation  (4.3.1)  to 

Cs  A’ 

H  - - ?s,  ’  (4.3.2) 

m 
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where 


( 4 , 3 .  J) 


Cs  =  t  ps  lae/apJ2. 

Equation  ( >i .  3 . 2 )  may  be  expressed  as  a  gradient  of  a 
potential  R*: 

?  R* ,  14.3.41 


where  R*  - - L  *  rs 

m  L 


5 

or  R  =  -  -  r  cosS 

m 


where  S  =  the  planet  centered  angle  between  the  Sun  and  the 
satellite.  This  angle  {see  occultation  subsection)  is 


cosS  =  CtcosQ  cosu  -  cosi  sinG  sinu) 

+  AtsinQ  cosu  +  cosi  cosQ  sinu) 

+  B(sini  sinu), 

where  A  =  cos6  sina 
B  =  sin6 
C  =  cos6  cosa 
6  =  declination  of  the  sun 
a  =  right  ascension  of  the  sun 
u  =  argument  of  latitude  (f+w). 

For  the  potential  above,  R  ,  in  the  classical  orbital 
element  set,  the  following  partial  derivatives  are  needed: 


(vector  form)  {  1 . 3 . 5 a  1 


(scalar  form)  (>1.3. 5b) 


Cs  A’  r 

-  —  cosS 

m  a 


^R* 


Cs  A' 


m 


£ r  ^(cosS) 

cosS  — —  +  r  — - - 

a  e  a  u 


<)  e 


(4.3.6a) 


(4.3.6b) 
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(4.3.6c) 


C>R* 

1l 

A*  ^(cosS) 

hi 

m 

C  *  J 

_  _  ^s_ 

A'  ^(cosSJ 

m 

b  Q 

^  R* 

-  _  El 

A*  $(cosS) 

^  Uf 

m 

*  W 

^  R*  =  Cs  A' 

b  H  b  f  b  M  m 

^ ( cosS )  - 

cosS  — —  +  r  — - 

ou  a  u 


( 4 . 3 . 6d ) 

( 4 . 3 . 6e  ) 

(4 . 3 . 6f ) 


These  further  partial  derivatives  are  needed  to  form  the 
soLar  radi.atj.on  equations  of  motion: 
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cosx  cosG  sinu) 


-cosS  =  Cl-cosQ  cosu  - 
3  G 

+  Al-cosQ  cosu  -  cosi  sinG  sinu)  (4.1.7 


=  Cl-cosG  cosu  -  cosi  sinG  cosu) 

5  u 

+  A( -sinG  sinu+cosi  cosG  cosul+B(sini  cosu),  (4.1.71 


where  the  partial  derivatives  below  have  been 
interchangeably.  (See  the  atmospheric  drag  section.) 

A,  =  JL  =  A. 

dw  3u 


The  new  orbital  equations  from  Equs.  (4.1.6)  are 


da  n2 

a 

0 

•# 

3c 

M  cosS ) 

dt  na 

r 

(i-e^J  cosS 

3  u 

f  r  ° 

d  u 

dh 

dt 


na^/i-e^ 
it/ l-e2 
na 


( sinu( 2+kcosu+hsinu )+h  )r 


cosu  cosS  - 


nk  coti 

2T  2 

na  i l-e 


^( cosS ) 
^  u 
^(cosS) 


di  =  ncoti  r  ^(cosS) 
dt  na2/l-e2  3  u 

dk  n 


^  i 

it  ^(cosS) 

na2  sin i  y’l-e2  &  Q 

^ (cosS ) 


dt  na2/l-e2 


(cosu( 2+kcosu+hsinu )+k)r 


3  u 


(4. 


(4. 


(4. 


) 

used 


4.8a) 


4.8b) 


4.8c) 
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iul-e2  rehcoti  r  ^(cosS) 

+  -  cosu  cosS  -  — - - —  — r - 

na  na  *  i-e“  c)  i 


l-l.-l.Bd) 


dfi  =  re  ^(cosSJ 

—  — -  c  sin i  ■— r - 

dt  na”< l-e  O  i 


(-l.-l.8eJ 


d  Ly  = 

nfJ  * 

I  0  A  l*  M  nil  A  U  C 

>  1  (-1  |  I  \  /  M  %  n  ||  __ 

^(cosS ) 

dt 

Syna2v l-e2 

\ itKcosutns 

anu )  \  Ksinu 

ncosu 1 r  * 

c>  u 

_  - f  a0‘/l-e2  { kcosu+hsinu)+2r 

S0na2L 

j  cosS 

» 

re 

Sq  -  cosi 

4(cosS)  +  n  d0g, 

_  /  1  1  0  r  \ 

T 

0  V 

Sgna^v l-e4 

sini 

—  -  -  \  “t  «  1  •  O  L  1 

SQ  dt 

where 

a 

11 

1 

0 

u 

a 

> 

1 

Occultation 

As  the  satellite  orbits  Mars,  it 

will  periodically 

swing 

behind  the 

planet  such  that  the 

Sun,  and  thus  its 

radiation,  is  blocked  or  occulted.  As  a  result,  we  must 
calculate  this  region  of  solar  occultation,  called  g(u). 

Consider  the  right  handed  coordinate  system  ( S *  Sg ,  S3  ) 
whose  origin  is  the  center  of  Mars.  The  axis  will  point 
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at  the  sun,  axis  S3  will  be  normal  to  the  plane  of  the 
equator,  and  S2  will  be  placed  so  as  to  compLete  the  right 
handed  system  as  in  Figure  4.1  below. 

tnder  this  system,  the  conditions  for  occultation  are: 

St  <  0  and  S22  +  S32  <  Rp2,  (4.4.91 

where  Rp=  radius  of  Mars  plus  atmospheric  altitude  blockage. 
Equation  (4.4.91  may  be  rewritten  using 

_2_c2.c2.q2 

r  -  +  S2  +  S3 

to  become 

r2  -  Sl2  <  Rp2  or  Sl  <  -  y'r2  -  Rp2.  (4.4.10) 


This  equation  requires  that  always  remain  more  negative 
as  the  condition  for  occultation. 


Figure  4.1:  Occultation  Geometry 
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Assume  the  orbital  elements  of  the  spacecraft  are  known 
and  the  right  ascension  (a)  and  the  declination  (6)  of  the 
viewed  object  are  known  in  the  same  coordinate  system.  The 
cartesian  position  of  the  spacecraft  in  standard  notation  is 
given  by 


=  r(cos£2  cosu  -  cosi  sinO  sinu) 
\2  =  r(sinQ  cosu  +  cosi  cosG  sinu) 
X3  =  r  sini  sinu, 


{•}.}.  lla) 
(4.4.11b) 


(4.4.  Ucl 


where  u  =  w  +  f.  In  order  to  express  the  cartesian 
coordinates  relative  to  the  -  coordinate  system,  two 
rotations  are  performed: 


cos{-6)  0  -sin(-6) 

0  1  0 
sin(-6)  0  cos{-6) 


cosa  sina  0 

-sina  cosa  0 

0  0  1 


Sj  cos(6)  cosa  cos{6)  sina  sina 

So  =  -sina  cosa  0 

S3  -sin(6)cosa  -sin(6)sina  cos(6) 


(4.4.12) 


This  is  now  used  in  equation  (4.4.11)  to  obtain 

Si  =  (cc  -sa)xi  +  (  cos6sina)x2+(  sin6 1x3  <-  v'r2  -  Rp2  (4.4.13) 

=  A  xL  +  B  x2  +  C  x3  <  -  yr2  -  R  2, 


where  the  variables  A,B,C  were  used  earlier.  Substituting 
for  the  cartesian  coordinates  changes  the  equation  to 


=  rCMcosG  cosu  -  cosi  sinG  sinu) 

+  rAtsinG  cosu  +  cosi  cosG  sinu) 

+  rBtsini  sinu)  £  -  <r^  -  Rp2  .  (>1.4.14) 

Redefine  the  coefficients  as  follows: 


F  =  CcosG  +  AsinG 
G  =  cosi(AcosG  -  CsinG)  +  Bsini 
H  =  vl-  (Rp/D2 


(4.  ). 15a) 


(4.4.15b) 


(4.4.15c) 


Finally,  substitution  makes  the  occultation  condition 


g(u)  =  Fcosu  +  Gsinu  +  H(u)  £  0. 


(4.4.16) 


The  occultation  condition  has  three  restrictions. 
First,  g(u)  depends  on  cosu  and  sinu,  causing  it  to  be 
periodic  with  period  2PI=2x3 . 14159 ... .  Next,  g(u)  also 
depends  on  H(u),  implying  the  minimum  of  g(u)  may  not  be 
zero.  Finally,  when  g(u)  is  at  a  minimum  with  u  =  um^n  and 
is  negative  the  satellite’s  orbit  is  occulted  at  least  part 
of  the  time  from  the  Sun.  Therefore,  to  check  occultation, 
one  must  find  the  minimum  of  g(u),  which  could  be 
accomplished  by  solving  for  the  roots  of  the  derivative  of 
g(u). 


57 


Itmaap  harts,  D.rag  Model 


The  following  model  applies  to  both  models  o  the* 
Martian  atmosphere  developed  earlier.  Tnoi>e  *?}del^ 
calculate  the  atmospheric  mass  density,  p,  to  i;-?  in 

conjunction  with  this  drag  model. 

The  force,  D,  per  unit  mass  due  to  atmospheric  drag  is 
given  by  (48:81): 

CdpAv2 

D  =  -  ,  (4.5.1) 

2m 

where  D  =  drag 
m  =  mass 

Cjj  =  drag  coefficient 
p  =  atmospheric  density 
A  =  effective  spacecraft  area  for  drag 
v  =  satellite  velocity. 

Using  the  notation  of  (30:81),  the  drag  acts  opposite 
to  the  velocity  of  the  satellite  and  may  be  written  in 
vector  notation  as 

12  =  -  D  up  (4.5.2) 

where  Up  is  the  unit  vector  along  the  tangent  in  the 
direction  of  satellite  motion. 

The  drag  force  is  broken  down  into  two  components: 
radial  and  transverse  components  of  the  satellite's 
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velocity.  These  are  defined  by  a  radial  unit  vector,  ur, 
along  the  position  vector  c  of  the  satellite,  and  the 
transverse  unit  vector,  uq,  perpendicular  to  r  and  forming 
an  acute  angle  with  the  velocity  vector  y  in  the  orbital 
plane.  Writing  the  drag  using  these  unit  vectors  yields 

C  =  Fr  ur  +  Fq  uq.  (>1.5.3} 

To  determine  these  drag  components,  Fr  and  Fq,  first 
define  the  velocity  vector  as 

Y  =  v  up  =  r  ur  +  rr  uq.  (4.o,4) 

This  makes  the  drag  components 

Fr  =  12  *  ur  =  -  D(Up  ‘  ur)  =  -  r  D  (4.5.5) 

v 

Fq  =  n  *  Uq  =  -  D(  Up  *  Uq  )  =  -  rjr  D  (4.5.6) 

v 

Making  use  of  these  standard  orbital  relationships  (12): 


a( 1-e2 ) 

r  =  - 

(4.5.7a) 

1+ecosf 

k  esinf  / M 
•  * 

r  =  .  .  9. 

(4.5.7b) 

/a( 1-e2 ) 


r2f  =  h  =  /k2Ma{  1-e2  ) 

(4.5.7c) 

k( 1+ecos  f)f M 

r'f  =  „ 

( 4 . 5 . 7d) 

v’a(  1-e2 ) 
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V 


I 


(  4 . 5 , 7  e  J 


kfM  i  ( l+e2  +  2ecosf) 

*a{ 1-e2 } 

by  substitution  into  equations  (4.5.5)  and  (4.1.6)  yields 
-esinf  D 

Fr  a  - t -  (4.5.8) 

i  (  l-!-eii+2ecos  f) 

-f l+ecos f )  D 

-  n  •  (4.5.9) 

/ ( l+e^+2ecos  f ) 

Now  place  these  magnitudes,  Fr  and  Fq,  into  the 
perifocal  orbital  reference  frame  specified  by  a  unit  vector 

A 

in  the  direction  of  perihelion,  P,  and  a  unit  vector  f=90* 

A 

in  the  direction  of  motion  from  P,  called  Q,  making  the  drag 

A  A  A  A 

H=  Frur+F0U0  =  (  Frcosf-F0sinf)P  +  ( Frsinf+F0cos  f  )Q  (4.5.101 
and  making  the  position  vector 

A  A 

r  =  (rcosf)P  +  (rsmf)Q 

a(l-e2)cosf  A  a(l-e^)  A 

=  -  P  +  -  Q.  (4.5.11) 

1+ecosf  ( l+ecosfjsinf 

The  last  step  is  to  calculate  the  partial  derivatives 
of  the  perturbation  due  to  atmospheric  drag  using  the 
following  equation: 

=  ia  ir  =  £  .  in  ,  (4.5.12) 

d  z^  dr  ^  z^ 
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where  the  z^  represents  the  six  classical  orbital  elements. 
The  resulting  partial  derivatives  are: 


=  Fr  r 


{  1 . 5 . 13a ) 


-  Fr  acos f  +  Fg  asinf  1+  - — 

a{  l-e“ ) 


( 4.5. 13b) 


^  R  =  Fr  ea  sinf  +  Fg  a2  yl-e* 
Am  t  i  ’  -e2 )  r 


(•1.5.13c) 


k  =  Fg  r  eosi 
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( -1.5.13d) 


=  Fg  r 


(4.5. 13e) 


(4.5. 13 f) 


Now  insert  these  into  Lagrange’s  Planetary  Equations  in 
order  to  find  the  orbital  equations  of  motion  in  the 
classical  orbital  element  set  that  are  due  to  drag: 


da  =  -  C^A  pna^  1+e  +  2ecosf 


(4.5.14a) 


2  "I  i  /o 

de  =  -  C(jA  pna  (e+cosf)  1+e  +  2ecosf  l/' 


d  w  =  - 


C^A  pna  sinf  [l+e  +  2ecosf 


(4.5.14b) 


(4.5.14c) 
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di  =  0 
dt 


(4.5.  14d) 


dll  =  0 
dt 


( 4 . 5 . 1 4  e  ) 


2  0 

dM  =  -  CjjA  una  sinf  (1+e  +6003/*)  vl+e *2ecos f 
dt  m  e  ( 1+ecos  f ) 


(4.5. 14 f) 


Convert  these  equations  into  the  modified  orbital  element 
set  and  substitute  to  get: 

*  *2  "  * 
dh  h  CjjA  una  {e+cosfl  1  +  e  +2ecosf 

dt  e  m  1-e^ 


-  C^A  una  sin f  1+e  +  2ecosf 


-  CjjA  una  (h+sinu)|l+e  +2ecosf  l1/2 


( 4 . 5 . 15a  J 


dk  _  k  -  C^A  pna  (e+cosf)  1+e  +  2ecosf 


dt  e 


-  CjjA  una  sinf  1+e  +  2ecosf 


=  -  C^A  una  ( k+cosu )  [*l+e^  +  2ecosf] 


(4.5.15b) 
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d£.j$  -  i  -  CfjA  pna  sinf*  (  l  +  e~+ecos  Hy  l+e22ecos f 
dt  Sn  m  e  (1+ecosf) 


C^A  pna  sinf  l+e"*+  2ecosf 


l  CjA  pna  esinf 

Sq  m 


P'  - 


y’l+e22ecos  f ,  (4.5.15c) 


ay'  l-e' 


where  ccosf=  kcosu  +  hsinu,  esinf=  ksinu  -  hcosu,  u =  v  +  f. 

Now  that  we  have  the  drag  equations  of  motion,  we  need 
to  apply  the  Method  of  Averaging  to  eliminate  fast  variables 
as  shown  in  reference  (28:3-26).  Retain  only  the  averaged 
effects  of  atmospheric  drag  over  one  period  in  order  to 
approximate  long  term  orbits.  The  averaged  rate  is  given  in 
(28:3-26)  as 


dX  = 

dt 


i  1 


p  d. £  dt, 
dt 


(*1.5.16) 


where  T  =  a,h,k,  L 

Time  is  related  to  the  true  anomaly,  f,  by 


dt  = 


r2  d  f 
na2/(l-e2) 


(4.5.17) 


Substitution  into  (4.5.16)  with  a  period  of  2PI  yields 


dr  n 
dt  2PI 


2PI  dr  r2df 
o  dt  na2/(l-e2) 


(4.5.18) 
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Over  one  period,  the  Argument  of  Periapsis,  w,  will  have 
only  a  small  change.  As  a  result,  we  can  write 

df  =  dv  +  du  si  du  or  df  -  du  (1.5.19) 


Substituting  the  orbital  equations  derived  above  into 
the  Integral  averaging  equation  yields: 

_  $  r  , 

da  =  2  [2PI  P  -C  ^  ( l+e2+2{  kcosu  +•  hsinu)]2/2du 

dt  l“e  JO  *  (4.5.20a) 


dh  =  *  f 
dt  J0 


x  ( l+e2+2{ kcosu  +  hsinu ) ) 1/,2du 


di  =  0 
dt 

dit  =  *t»  f2PI  p  r  2  ( k+cosu ) 

dt  J  o  « 

x  t  l+e2+2{  kcosu  +  hsinu )]  ^^du 


da  =  0 

dt 


(4.5.20b) 

(4.5.20c) 

(4.5. 20d ) 
(4.5. 20e ) 


di«M  =  *  f2PI  p  £ 
dt  j0  a 


2  ( ksinu-hcosu ) [ l+e2+2( kcosu+  hsinu))1/2 


r 

a / ( 1-e2  ) 


(4.5. 20f ) 
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where  «t>  = 


L 


2 

-  C^A  na 

m  1-e2  2PI 

and  h  is  calculated  using  one  of  the  two  Martian  atmosphere 
models  that  were  explained  in  section  III. 


Third  Body  Disturbing  Potential  (Solar  Gravity} 


The  derivation  of  these  equations  is  quite  extensive. 
Further  information  may  be  found  in  reference  (28)  and  other 
works.  The  following  equations  are  those  necessary  for  this 
model . 


The  averaged  third  body  disturbing  potential  R*  is 
given  as  follows: 


p 1  • 

R'  =  —  S 


a  ' 


n=2 


n  v  v  9  fp 

m=0  Kn®  p=oL 


nmp(^  ^npq^e^ 


] 


x  hI0  Fnmh^’ J  .  ^[Gnhjte’ )  cos'pjj  , 


(4.6.1) 


where  o  =  -qw  -  (n-2h)w'  -  (n-2h+j)M’  +  mtG-Q’),  (4.6.2) 

^nu  “  ii-LU~a)  !.  i  KnQ  =  1,  q  =  2p  -  n , 

(n+«) ! 


H  is  defined  latert  and  the  variables  with  the  superscript  * 
refer  to  those  of  the  disturbing  body. 

Rewriting  the  disturbing  potential  R’  using  definitions 
from  (28:3-21)  obtains: 
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R'  =  a,  “  m=0  Knt"  pio[FnraPti}  H'npq(e) 

H“2 


a  p 

x  K?n  ^nmh^  *  “  Gnhj*e  ^  ^nmqhj 

j=-«  L 


(1.6.3) 


The  function  Hnp(2p-n)  *s  defined  as 


‘np( 2p-n )  = 


( -p }n“2p  * 
(  i+02)n+1 


2n+i-2p‘  S  n+1  2p'+l 

n-2p>  E  k  it 


(4.6.4) 


n-2p ' +k 
k 


where  p'  =  p,  p  <  n/2 

P'  =  n-p,  p  >  n/2 

k  =  all  integers  where  binomial  coefficients  are 
nonzero . 


Therefore,  the  partial  derivatives  of  R*  in  the 
modified  orbital  element  set  are: 


jM*’  =  U’  1  S  n  n  mSQ  Knn  |Q  Fnrap(i)  H'npq(e) 
da  a  a  n— 2  ® 


*  h=0  Rnmh^’^.  ^  Gnhj^e’*  Gnmqhj 
J=“® 


(4.6.5a) 
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<9 


=  u 

a'  n=2 


n 


n  .  n 

m=0  K°m  p=0  Fnmp{i) 


h  dH’npq(e 
e  de 


x  h^O  Fnmh^  *  “  ^nhj*e'^  ^nmqhj 

J =~" 

+  "’npql'l  hl0  Fnmhti’)j °nhjl<=‘>  ^ 


(  I  . 


&'  =  u'  s 

di  a*  n=2 


n 


v  K 
-«  >» 


n  ^nmp*  *■  ) 


m=0  ‘ nn  p=0  di 


^ 1 npq*  e  * 


X  h=0  Fnmh^  ^  Gnh j ^ e ' ^  Gnraqhj 
.i=-« 


(4. 


is1  =  U’  S 


d  k 


n=2 


“ 

1 

•sT* 
* _ 

n 

m=0  K°ra  p=0  Fnmp<l) 

i  j 

-  dH’*We 

e  de 


X  h=0  Fnmh*  Gnh j ^ e  '  ^  Gnmqhj 

J 


+  H’npqfe*  hg0  Fnnh‘l*»  2  W>'1  ~T 

J  =-*  O 


d  snmqhj 


(4. 


•  r 


i  - 


~  U*  2 
n=2 


ra?0  Knm  pIo  Fnmp^i)  H’npq(e} 


n 


X  h=0  ^nmh^ 1  ^  ^ 

j=-® 


Gnhj(  =  ') 


(4 


6.5b  J 


6.5c) 


) 


6 . 5d ) 


6 . 5e ) 
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iB'  =  0. 


( 4 .  6 . 5  n 


All  that  remains  to  be  done  is  to  substitute  into  the 
Lagrange's  Planetary  Equations  (-1.1.4)  as  we  have  done  with 
the  other  disturbing  potential  functions,  replacing  R  with 
R'  to  form  the  differential  equations  for  the  third  body, 
the  Sun.  These  substitutions  will  not  be  done  here. 
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ill.  Approach  and  Results 


Apg.r.oacii 

\  variety  of  assumptions  are  made  in  constructing  the 
dynamic  and  static  atmosphere  models.  Host  of  these  have 
already  been  explained  in  the  theory  sections.  This  section 
will  review  the  inputs  to  the  computational  models  based  on 
the  theoretical  models  of  the  previous  sections. 

The  Martian  geopotential  was  modeled  using  only  Jg  and 
,  the  dominant  terms  representing  the  strong,  first-order 
equatorial  oblateness  effect  and  significant  first-order 
asymmetry  of  Mars  (28:2-35  )  and  (8:1).  (See  Table  5.1.) 
[n  addition,  the  orbital  missions  were  chosen  close  to  a 
known  solar  epoch  and  Prime  Meridian  (28:2-10)  in  order  to 
minimize  error  due  to  the  coordinate  system  used,  namely  the 
planet  equator  of  epoch  coordinate  system  for  Mars. 
Finally,  empirical  data  was  used  to  choose  limits  for  the 
atmosphere,  both  the  "outer"  limit  and  the  "inner"  limit. 

The  "outer"  limit  is  set  to  1000  km  since  this  is  where 
the  first  "heavy"  ions,  atomic  oxygen,  exist  in  significant 
concentrations  (37:210).  Although  atomic  hydrogen  has  been 
observed  as  high  as  10,000  km,  these  are  much  "lighter" 
atoms  and  so  exert  far  less  drag.  Hence,  the  1000  km  limit. 
The  "inner"  limit  is  set  to  130  km  based  on  the  fact  the 
ionosphere  constituent  concentration  is  greatest  at  this 
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altitude  (21:27).  In  addition,  the  turbopause  altitude  (the 
dividing  line  between  the  denser  lower  atmosphere  and  the 
less  dense  ionosphere)  is  very  close  to  this  region.  It 
may  be  expected  that  any  satellite  which  reaches  this 
altitude  is  actually  beginning  to  enter  the  atmosphere 
proper  and  will  shortly  impact  Mars  due  to  t  lie  much 
thicKer  atmosphere. 

The  mean  anomaly  M,  which  is  a  fast  variable  (short 
term  effects)  is  averaged  out  in  the  theoretical  models.  As 
a  result,  it  can  assume  any  value.  For  convenience,  it  has 
been  chosen  to  be  zero  and  kept  constant.  The  other  orbital 
parameters  (a,e,i,2,w)  could  be  varied  at  will.  The 
initial  values  were  based  on  data  from  the  Viking  missions 
(■11:3979-3982),  Mariner  9  (26:484  ),  and  from  scientific 
missions  to  Mars  (45)  and  (28).  A  representative  sampling 
of  data  for  possible  real  life  missions  was  obtained, 
covering  both  close  and  far  orbits  as  well  as  both  near 
equatorial  and  near  polar  orbits. 

The  effects  on  orbital  motion  of  using  a  dynamic 
Martian  atmosphere  model  are  demonstrated  by  comparing  the 
resulting  orbital  element  data  from  the  dynamic  model  to 
those  from  the  static  model.  The  orbital  parameters 
examined  were  a,  e,  w,  and  the  periapsis  altitude,  hp.  The 
periapsis  altitude  is  calculated  as  follows: 

hp  =  a(l-e2)-rm,  (5.1) 
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where  rm  is  the  reference  altitude  of  the  surface  of  Mars, 
which  is  calculated  differently  by  the  two  atmosphere 
models.  Both  sets  of  data  are  graphed  as  a  function  of 
time.  The  Argument  of  Periapsis  is  only  graphed  for  the 
first  four  cases  in  order  to  check  for  consistency  with  the 
ocher  orbital  elements  and  for  any  trends. 

The  satellite  model  requires  input  values  for  the 
satellite  mass,  the  satellite  drag  coefficient,  and  the 
effective  satellite  area,  which  is  needed  for  both  drag  and 
solar  radiation  calculations.  Typical  values  were  chosen 
for  these  variables  from  references  ( 23 ) , ( 28 ) , ( 41 ) , and  (45). 
Table  5.2  lists  these  values. 

The  theoretical  models  were  programmed  using  the 
FORTRAN  programming  language  and  combined  programs  both 
self-written  and  imported  from  (28)  and  (8).  Each  computer 
integration  run  used  a  step  size  of  5  days,  unless  otherwise 
noted,  and  were  typically  set  for  a  two  year  run  as  that  was 
deemed  sufficient  to  show  the  desired  effects  on  the  orbit. 
Special  attention  was  paid  to  low  orbits  to  verify  the  130 
km  altitude  limit  of  the  "real"  Martian  atmosphere. 


Table  5.1:Martian  Geopotential 

Coefficients  for  J2  and  J3 

n 

m  1  cnm 

1  ^nm 

2 

0  |  -0. 1960387250000000E-02 

|  0.0000000000000000E+00 

3 

|  0  |  -0.30634 19400000000 E- 04 

j  0.0000000000000000E+00 
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Table  5.2:  Spacecraft  Data  Input  Values 


Mass  (kg) 

i 

1000 

Drag  Coefficient 

O 

CM 

2 

Area  for  Drag  (km  ) 

10.0E-6 

Area  for  Solar  Radiation 

0 

Pressure  ( km“ )  : 

20.0E-6 

Solar  Radiation  Pressure 

o 

Constants  kg/kmsec~ ) : 

1.95E-3 

Table  5.3:  Other  Necessary  Input  Data 

Prime  Meridion  ( *  ) 

333.55971 

Solar  Epoch 

19910524  =  24 

May  1991 

Initial  Propagation  Time: 

19911007  =  7 

October  1991 

GMmars  Um3/sec2) 

4 . 2828287E4 

Radius  of  Mars  (km) 

3397.2 

Rotation  of  Mars  {‘/sec) 

4 . 061249803E- 

3 

Radius  +  Atmospheric  Blockage: 

Radius  +  90  = 

3487.2  (km) 

GMjun  (km3/sec2) 

0. 13271244E12 

Each  data  run  involved  73  points  for  each  parameter  for 
each  year  of  the  data  run.  Hence (  it  was  not  practical  to 
show  every  point  in  tabular  form,.  The  following  data  tables 
show  the  initial  parameters  at  the  time  equals  zerp  row  and 
then  a  data  point  for  each  120  days  of  the  run  for  as  long' 
as  the  run  lasted.  The  first  two  cases  are  shown  in  this 
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chapter  in  both  tabular  and  graphical  form;  the  other 
fifteen  cases  are  presented  in  Appendices  A  and  B. 

Results 

The  data  indicate  there  is  a  difference  in  the  orbital 
motion  of  a  satellite  due  to  a  dynamic  Martian  atmosphere  as 
compared  to  motion  resulting  from  a  static  Martian 
atmosphere.  The  dynamic  atmosphere  provides  lesser  drag  on 
the  spacecraft,  implying  a  lesser  atmospheric  mass  density, 
i.e.,  motion  through  the  static  atmosphere  tends  to  spiral 
the  satellite  into  the  surface  of  the  planet  more  quickly 
than  the  dynamic  atmosphere. 

This  effect  is  most  evident  for  low  orbits —  orbits 
only  a  few  hundred  kilometers  above  the  surface  of  Mars.  As 
the  orbital  ellipse  increases  in  size  (as  measured  by  the 
length  of  the  semimajor  axis),  the  amount  of  differentiation 
between  the  two  models  decreases  in  magnitude.  For  large 
orbits,  those  on  the  order  of  tens  of  thousands  of 
kilometers  in  the  length  of  the  semimajor  axis,  the 
differences  are  negligible.  (See  cases  8,10,11,3,  for 
example.)  Furthermore,  the  data  indicate  the  variation 
between  the  two  models  increases  with  time.  (See  cases 
1, 4,1*1.)  The  larger  the  size  of  the  orbit  the  longer  the 
two  models  remain  essentially  the  same  before  beginning  to 
differentiate.  (See  cases  6,9,10,8.)  Nonetheless,  once  the 
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differentiation  begins,  it  increases  as  time  goes  on.  In 
other  words,  the  amount  of  variation  between  the  two  models 
depends  on  both  orbital  size,  as  measured  by  length  of  the 
semimajor  .in is,  and  on  time.  This  is  logical  since  the 
larger  the  orbit  the  less  time  the  spacecraft  spends  m  the 
atmosphere,  resulting  In  less  change  in  the  orbit  due  to 
either  model,  and  hence  the  lesser  degree  of  variation 
between  the  two.  Once  the  two  trajectories  have  begun  to 
vary,  the  orbit  with  the  greater  atmospheric  drag,  always 
the  static  case,  will  cause  the  spacecraft  to  spend  a 
marginally  longer  amount  of  time  in  the  atmosphere,  thereby 
increasing  the  amount  of  variation. 

Varying  the  inclination  produces  no  apparent  effect  on 
the  orbits.  (See  cases  9,14,17.)  Furthermore,  changing  the 
eccentricity  while  keeping  the  other  orbital  elements  as 
they  were  only  has  an  impact  in  that  higher  eccentricities 
reduce  the  periapsis  altitude  and  hence  the  satellite  spends 
a  greater  amount  of  time  in  the  atmosphere.  (See  cases  7,8.) 
As  has  already  been  mentioned,  the  amount  of  time  spent  in 
the  atmosphere  is  the  crucial  determinant  in  variation 
between  the  two  models.  In  addition,  the  various  orbital 
parameters  agree  internally  with  each  other.  For  example, 
the  eccentricity  for  the  static  model  will  become  more 
circular  faster  than  the  eccentricity  of  the  dynamic  case, 
thus  indicating  the  greater  amount  of  drag  affecting  the 
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satellite  In  the  static  case  and  agreeing  with  the  decrease 
in  periapsis  altitude  and  seraimajor  axis  length. 

The  expected  beginning  of  the  "real"  atmosphere,  called 
the  lower  limit  previously,  was  at  about  130  km.  This  has 
been  verified.  For  those  trajectories  that  "blow  up,"  the 
data  show  that  once  the  periapsis  altitude  falls  below 
200  km  the  slope  of  the  curve  increases  to  nearly  vertical. 
The  trajectory  then  continues  until  about  130  km  ±  15  km 
where  the  numerical  calculation  "blows  up,"  thus 
indicating  that  the  satellite  impacted  the  surface.  (See 
cases  2,4.)  Doth  models  fall  within  this  range.  However, 
the  static  model  reaches  the  "blow  up"  or  impact  point 
first.  (See  case  5,  for  example.) 

In  summary,  the  dynamic  model  exerts  less  drag  on  the 
spacecraft  than  the  static  model.  The  variation  between  the 
two  models  increases  with  both  time  and  decreasing  orbit 
size.  Also,  the  orbital  parameters  agree  with  each  other. 
Finally,  the  limit  of  the  "real"  atmosphere  or  lower  limit 
is  shown  to  be  in  the  expected  region,  about  130  km. 


However,  a  closer  look 

should 

be  taken 

at 

the 

difference  in  the  atmospheric 

density 

between 

the 

two 

models.  This  difference  is  plotted  in  Fig.  B.48  for  a  two 
year  Case  1  input.  As  can  be  seen,  the  dynamic  density 
varies  in  an  oscillatory  manner,  which  is  expected  from  the 
confused  and  variable  state  of  the  Martian  atmosphere.  The 
static  density  varies  smoothly,  which  is  expected  due  to  the 
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exponential  function.  However,  while  the  dynamic  density  is 
sometimes  greater  than  the  static  density,  more  often  it  is 
less,  particularly  in  the  latter  half  of  the  run.  The 
static  case  was  modeled  using  a  reference  density  that  was 
an  upper  bound,  6  jc  10  kg/m  .  For  comparison  purposes, 
choose  another  typical  value  from  the  empirical  data  for  the 
height  region  of  the  orbits;  keeping  the  same  order  of 
magnitude  but  reducing  the  multiplying  factor  from  6  to  3.1 
modifies  the  relationship  between  the  two  models  (Figures 
B.-19  to  B.52).  As  can  be  seen,  the  two  models  are  now  much 
closer  to  each  other.  Nonetheless,  over  a  long  term  period, 
for  this  case  exceeding  one  year’s  run,  the  basic  difference 
persists:  the  static  case  will  impact  Mars  first  since  the 
static  density  is  still  higher  over  the  longer  time  period 
(see  Fig.  B.52).  However,  if  the  static  reference  density 
is  further  decreased,  the  lifetime  of  the  static  satellite 
will  continue  to  improve  relative  to  the  dynamic  satellite. 
But,  until  empirical  data  provides  a  better  agreed-upon 
value  for  the  reference  density,  the  best  that  can  be  said 
is  that  the  dynamic  model  exerts  less  drag  over  long  term 
orbits. 
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Table 

5.4:  Case  1 

Data 

Dynamic 

Data 

(5  year  run) 

t 

a(  kra ) 

e 

vi  *  ) 

hp ( km) 

0 

3758.31700 

.00704800 

270.00000 

334.62838 

3757.92726 

.00704882 

269.90415 

334.23831 

210 

3756 . 89135 

.00706625 

269.85892 

333.14423 

360 

2755 . 90921 

.00707127 

269.93376 

332.15018 

•180 

3755 . 6384  1 

.00705635 

269.99388 

331.93731 

600 

3755 . 52698 

.00704709 

269.92522 

331.86143 

720 

3755.42520 

.00705840 

269.85957 

331.71790 

840 

3755 . 16883 

.00707317 

269.93145 

331.40790 

960 

3754 . 60034 

.00706864 

270.00108 

330.86043 

1080 

3754 . 20377 

.00705504 

269.93150 

330. 51773 

1200 

3754 . 12023 

.00705701 

269.85583 

330.42736 

1380 

3754.06216 

.00706860 

269.96377 

330.32618 

1500 

3753. 96388 

.00706271 

270.01380 

330.25074 

1620 

3753. 66256 

.00706013 

269.92361 

329.9S123 

1740 

3753.39308 

.00706717 

269.84392 

329.66720 

1800 

3753.35281 

.00707068 

269.86774 

329.61405 

Case  1 

Continued  (Static  Data) 

0 

3758.31700 

.00704800 

270.00000 

334.62838 

120 

3756.87252 

.00706130' 

269.98956 

333.14409 

240 

3755.36809 

.00707024 

270.03615 

331.61674 

360 

3753.79840 

.00706500 

270.04730 

330.07780 

480 

3752.15764 

.00706019 

270.01820 

328.46670 

600 


3750.43952 


00706589 


270.00703  326.73934 


840 

3746.73809 

.00708106 

270.08833 

323.00720 

960 

3744.73504 

.00708230 

270.06910 

321.01371 

1080 

3742.61448 

.00708324 

270,02668 

318.90464 

1200 

3740.36238 

.00708801 

270.04725 

316 . 55067 

1320 

3737 .96176 

.00709282 

270.10459 

314.24908 

1440 

3735.39101 

.00709669 

270.13623 

311.68210 

1560 

3732.62321 

.00710448 

270.12866 

308.90486 

1680 

3729.62608 

.00711197 

270.09812 

305.90109 

1800 

3726.35802 

.00711701 

270.11560 

302.63749 

Table  5.5:  Case  2  Data 


Dynamic  Data 


3697.20000  .00250000 


w(  ‘  )  hp(  km) 


270.00000  290.75700 


120 

3695.72641 

.00528445 

229.39501 

278.99653 

240 

3690.66345 

.00876479 

237.82332 

261.11556 

360 

3684.48252 

.01097531 

252.23568 

246.84420 

480 

3682.96559 

.01184056 

267 . 35057 

242.15722 

600 

3682.27339 

.01142396 

282.36780 

243.00723 

720 

3681.66220 

.00982776 

296.37131 

248.27977 

Case  2 

Continued  (Static  Data) 

0 

3697.20000 

.00250000 

270.00000 

290.75700 

120 

3687.99932 

.00518093 

236.25771 

271.69205 

240 

3675.91680 

.00763305 

243.72583 

250 .65834 

360 

3658 .04837 

.00862316 

254.07183 

229.30444 

480 


3622.10189 


00824591 


263.24949  195.03439 
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5.5:  Periapsis  Altitude  vs.  Time  (Case 
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Time  (Days) 


200.00 


VJLl  Conclusions  and  Recommendations 


Conclusions 

The  effect  on  long  term  orbital  motion  about  the  planet 
Mars  using  an  initial  dynamic  model  of  the  Martian 
atmosphere  was  studied.  The  effect  is  noted  through 
comparing  the  orbital  elements  generated  by  the  dynamic 
model  with  those  generated  by  a  static  model  of  the  Martian 
atmosphere.  The  following  conclusions  are  evident: 

L.  The  dynamic  model  exerts  less  atmospheric  drag  on  the 
satellite,  implying  a  less  dense  atmosphere  and  increasing 
the  lifetime  of  the  satellite. 

S i.  The  variation  between  the  models  increases  with 
decreasing  orbit  size,  as  measured  by  semimajor  axis  length 
for  fixed  periapsis  altitude,  and  length  of  time  of  the 
trajectory.  A  further  factor  is  the  input  data  chosen  for 
the  static  model,  particularly  the  reference  density;  the 
smaller  the  reference  density  the  smaller  the  amount  of 
variation  between  the  two  models. 

j.  Variation  of  other  orbital  parameters  than  the  semimajor 
axis  shows  only  limited  effect  on  the  orbit  unless  periapsis 
altitude  is  decreased  as  a  result  of  the  change. 
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!•  The  orbital  parameters  were  internally  consistent,  which 
indicated  the  model  worked  correctly  and  supported  the 
conclusions  across  severaL  parameters. 

5.  The  location  of  the  lower  limit  of  the  atmosphere,  the 
beginning  of  the  dense  part  of  the  atmosphere,  was  verified 
to  be  at  approximately  130  km  altitude. 

KecQnuaen.dac.iQU3 

Based  on  the  results  of  this  investigation  and  on  the 
research  that  was  used  to  form  the  models,  the  following 
recommendations  for  further  study  are  proposed: 

1.  Further  data  analysis  of  the  space  probe  data  from  Mars 
to  develop  a  better  understanding  of  the  Martian  atmosphere 
in  order  to  facilitate  model  development  of  the  various 
effects  present  in  the  atmosphere  as  well  as  to  provide 
standard  numbers  for  such  inputs  as  the  reference  density  in 
the  static  model. 

2.  Development  of  a  model  which  includes  a  greater  number 
of  the  effects  existing  in  the  Martian  atmosphere;  e.g., 
polar  cap  sublimation/condensation,  atmosphere  wave 
perturbations,  and  latitude/longitude  gradients. 
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«J.  Further  research  generating  numerical  data  using  this 
and  subsequent  models  along  with  an  in  depth  analysis  of  the 
resulting  data  in  order  to  gain  a  large  data  base  for  the 
detailed  evaluation  of  dynamic  Martian  atmosphere  effects  on 
long  term  orbits. 

•1.  Form  an  atmosphere/trajectory  model  that  includes  short 
term  orbits  and  effects  in  order  to  more  realistically  and 
accurately  represent  the  Martian  atmosphere's  effect  on 
orbits . 
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Appendix  A 


Supplemental  l&blca  far  Section*  1  *nd  £1 

The  following  pages  contain  tabulations  of  the  data 
generated  for  cases  three  through  seventeen. 


90 


Table 

A.l:  Case 

3  Data 

Dynamic 

Data 

1 

a  ( km ) 

e 

w(  } 

hp( km J 

0 

3725.00000 

.00700000 

35,00000 

301.72500 

120 

3721 . 12788 

.01383227 

311.92456 

275.41472 

■n 

3721.07098 

.01524739 

252 .85715 

267.13436 

360 

3717.80320 

.01010305 

190,60954 

283.04203 

•180 

3716.90853 

.00538407 

73 . 53509 

299.69644 

600 

3716.61497 

.01130528 

334.78240 

277.39759 

720 

3716.32255 

.01512161 

277.59429 

262.92579 

Case  3 

Continued  (Static  Data) 

0 

3725.00000 

.00700000 

35.00000 

301.72500 

120 

3719.66732 

.01308124 

312.33641 

273.80945 

240 

3713.78845 

.01407112 

258.51358 

264.33127 

360 

3706.99647 

.00961410 

209.89607 

274.15703 

480 

3698.84002 

.00323519 

145.66006 

289.67357 

600 

3688.71651 

.00384913 

359.54218 

277.31814 

720 

3675.32101 

.00694169 

318.15801 

252.60807 

Table  A. 

2:  Case  4 

Data 

t 

a(  km) 

e 

v(*  ) 

hp{ km) 

Dynamic 

Data 

0 

3700.00000 

.00700000 

35.00000 

276.90000 

120 

3698.12682 

.01223392 

330.14531 

255.68423 

240 

3692.44248 

.01532293 

286.81561 

238.66345 

360 


3682.67101 


01516099 


251.64807  229.63806 


480 

3680.42219 

.01278364 

219.20308 

236.17300 

600 

3679.60498 

.00906081 

181.90217 

249.06478 

720 

3678.71460 

.  00553958 

124 .75665 

261.13606 

Case  4 

Continued  (Static  Data} 

0 

3700.00000 

.00700000 

35.00000 

276.90000 

120 

3688.67503 

.01064301 

330.23249 

252.21642 

240 

3673.72708 

.01183673 

295.34194 

233.04217 

360 

3650.31611 

.01051592 

276.89075 

214.72966 

480 

3581.29682 

.00692021 

272.61030 

159.31348 

Table 

A. 3:  Case  5 

Data 

Dynamic 

Data 

t 

a(  km ) 

e 

w(*  ) 

hp(km) 

0 

3650.00000 

.00700000 

35.00000 

227.25000 

120 

3639.12449 

.00770013 

21.16953 

213.90275 

200 

3534.94210 

.00583657 

49.89566 

117.11016 

Case  5 

Continued  (Static  Data) 

0 

3650.00000 

.00700000 

35.00000 

227.25000 

120 

3537.30945 

.00219324 

251.91097 

132.35129 

Table  A. 4:  Case  6 

Data 

Dynamic 

Data 

t 

a(km)  e 

*(  *  ) 

hp(km) 

0 

3700.00000  .00700000 

35.00000 

276.90000 

15 

3699.90603  .01225080 

209.15020 

257.37920 

25 

3699.69948 

.01580484 

269.62480 

244.02631 

1600 

3606.31746 

.01280540 

309.01570 

162.937)4 

Case  6 

continued  (Static  Data) 

0 

3600.00000 

.00700000 

35.00000 

...  .i  .  < J-  • 

177.60000 

15 

3578 .35719 

.01131349 

220.85973 

140.67348 

25 

3542 .95578 

.00947071 

276.84234 

1 12.20148 

Table  A. 

5:  Case  7 

Data 

Dynamic 

Data 

t 

a(  km ) 

e 

v(  ) 

hp( km) 

0 

40850.00000 

.89730000 

35.00000 

798.09500 

120 

40850.17026 

.90239985 

33.76714 

589.78284 

240 

40850.15743 

.90189507 

33.40832 

610.40202 

360 

40850 . 08705 

.90395227 

32.37761 

526 . 35814 

480 

40849.87620 

.90793429 

31.92682 

36' .67302 

600 

40849.86727 

.90561511 

31.69218 

458.41004 

720 

40850.04342 

.90774866- 

31.05805 

371.27116 

Case  7 

Continued  (Static  Data) 

0 

40850.00000 

.89730000 

35.0000 

798 . 09500 

120 

40850.17025 

.90239985 

33.76714 

589.78282 

240 

40850.15485 

.90189506 

33.40832 

610.40200 

360 

40850.08699 

.90395230 

32.37762 

526.35672 

480 

40849.66663 

.90793384 

31.92682 

363.67195 

600 

40849.46672 

.90561426 

31.69216 

458.40713 

720  40849.57030 


90774763 


31.05802  371.26966 


Table  A. 6:  Case  8  Data 


Dynamic 

Data 

mm 

a(  km) 

e 

w(  ) 

hp( km ) 

0 

•10850.00000 

.88010000 

35 . 00000 

1500.71500 

120 

■10850.14119 

.88547510 

34.06634 

1281.15819 

2-10 

•10850 . 14505 

.88513472 

33.97313 

1295.06333 

360 

40850.14171 

.88688502 

33.11451 

1223.56290 

■180 

40849.93049 

.89146582 

32.82584 

1036.41370 

600 

40849.92690 

.88908986 

32.69441 

1133.47097 

720 

40850.03299 

.89062490 

32.09720 

1070.77649 

Case  8 

Continued  (Static  Data) 

0 

40850.00000 

.88010000 

35.00000 

1500 . 71500 

120 

40850.14149 

.  8854  7510 

34.06634 

1281.15819 

2-10 

40850 . 14503 

.88513472 

33.97313 

1295.06333 

360 

40850.14171 

.88688502 

33.11451 

1223.56290 

■180 

40849.93049 

.89146582 

32.82584 

1036.41370 

600 

40849.92690 

.88908986 

32.69441 

1133.47097 

720 

40850.03299 

.89062490 

32.09720 

1070.77649 

Table  A 

l.7:  Case  9 

Data 

Dynamic  Data 

HU 

a(  km) 

e 

v(*  ) 

hp( km) 

0 

20000.00000 

.81500000 

35.00000 

302.80000 

120 

20000.14489 

.81707038 

30.55843 

261.41888 

mat 

20000.09376 

.81652226 

26.38588 

272.37194 

360 

19999.88509 

.81900112 

22.13464 

222.75686 

■180 

19999,88517 

.81914910 

18.15247 

219.79727 

MJU 

19999.99772 

.81970941 

14.17258 

208.61129 

720 

20000.10522 

.82134796 

10.25674 

175 . 85950 

Case  9 

Continued  (Static  Data) 

0 

20000.00000 

.81500000 

35.00000 

302.80000 

12U 

19997.16488 

.81704281 

30.55800 

261.42503 

240 

19991.53239 

.81644413 

26.38310 

272.36303 

360 

19983.11872 

.81884697 

22.12759 

222.80246 

480 

19962.07634 

.81880413 

18.13516 

219.84581 

600 

19942.69444 

.81918958 

14.13833 

208.64696 

720 

19896.59695 

.82040012 

10.19529 

176.22649 

Table  A. 8:  Case  10  Data 


Dynamic  Data 


a(  km ) 


v(  ‘  ) 


hp { km ) 


0 

10000.00000 

.63530000 

35.00000 

249.80000 

120 

10000.01059 

.63637648 

19.28367 

239.03901 

240 

9999.72140 

.63753083 

3.64870 

227.39073 

360 

9999.63490 

.63849873 

348.09572 

217.70236 

480 

9999,59779 

.63887910 

332.57364 

213.8637 9 

600 

9999.49199 

.63926847 

317.06785 

209.93206 

720 

9999.68040 

.63922917 

301.51692 

210.39302 

Case  10  Continued  (Static  Data) 

0  10000.00000 


63530000 


35.00000  249.80000 


120 


9993.83322 


.63615189 


19.27510 


239.03732 


2-10 

9985.15170 

.63700294 

3.61037 

227.38074 

360 

9972.37600 

.63750990 

347.99307 

217 . 68755 

-180 

9957.81858 

.63736709 

332.36542 

213.83274 

600 

9943.16039 

.  63723235 

316.70995 

209.85690 

720 

9929.67438 

.63671780 

300.96585 

210.07400 

Table  A. 9:  Case  11  Data 


Dynamic  Data 


a  (km) 


5000.00000  .27060000 


K-C  ) 


35.00000 


hp( km) 


2-19.80000 


120 

4999.82315 

.27639218 

321.55416 

220.71115 

240 

4999.24985 

.27789333 

249.28687 

212.79168 

360 

4995.78748 

.27247961 

176.40634 

237.33727 

480 

4995.71358 

.26798495 

101.86486 

259 . 73753 

600 

4995.70813 

. 27050704 

26.99135 

247.13392 

720 

4995.70066 

.27625108 

313.53795 

218.43295 

Case 

11  Continued 

[Static  Data) 

0 

5000.00000 

.27060000 

35.00000 

249.80000 

120 

4993.73020 

.27552065 

321.45260 

220.65443 

240 

4985.17939 

.  27584452 

248.80682 

212.84497 

360 

4977.42737 

.26963612 

175.28874 

238.13316 

480 

4973.82942 

.26467018 

99.90200 

260.20509 

800 

4971.05403 

.26705131 

24.12279 

246.32757 

720 

4964 .24556 

.27183542 

309.66514 

217.58779 

Table  A. 10:  Case  12  Data 


Dynamic  Data 


:  a (km) 


5000.00000  .27060000 


hp( km) 


35.00000  219.80000 


120 

4999.60428 

.27082420 

13.94098 

248.39044 

240 

1999.07370 

.27113772 

353.17342 

246.43624 

360 

4996.34753 

.27109197 

333.39445 

244.67784 

480 

4996.14283 

.27130026 

314.36703 

243.48798 

600 

4996.10494 

.27144986 

295.41280 

242.71297 

720 

4996.03811 

.27155605 

276.47089 

242.13372 

Case 

12  Continued 

(Static  Data) 

1  1 

5000.00000 

.27060000 

35.00000 

249.80000 

120 

4995.35223 

.27018277 

15.18149 

243.49412 

240 

4990.66638 

.26982746 

357.96249. 

246.84755 

360 

4985.98804 

.26937066 

343.36737 

245.70918 

480 

4981.33806 

.26885673 

331.41692 

244.87178 

600 

4976.67017 

.26829597 

322.13487 

244.24964 

720 

4971. 96437 

.26771465 

315.53546 

243.69665 

Table  A. 11: 

Case  13  Data 

Dynamic  Data 

■ 

v(*  ) 


hp( km) 


0 

10000.00000 

.63530000 

35.00000 

249.80000 

120 

10000.15440 

.63498769 

259.90822 

252.97946 

240 

10000.21016 

.63393996 

124.93690 

263.47730 

360 

10000.10408 

.63521187 

351.60167 

250 , 71925 

180 

9999.74419 

.63609138 

216.29255 

241.79315 

600 

9999.60265 

.63494866 

83.38914 

253.16839 

730 

9999.60450 

.63581136 

309.91811 

244.51236 

Case  13 

Continued  (Static  Data) 

0 

10000.00000 

.63530000 

35.00000 

2  19.80000 

120 

9994.95423 

.63479912 

260.02277 

252.96605 

210 

9990.65152 

.63359845 

125 .37798 

263.39020 

360 

9985.43387 

.63468131 

352 . 58220 

250 . 66565 

*180 

9979.11712 

.63534572 

218.06626 

241.72776 

600 

9972.40889 

.63394591 

86.16906 

253.24105 

720 

9966.37191 

.63460121 

314.11492 

244.50021 

Table  A. 12 

:  Case  14  Data 

Dynamic  Data 

BS 

a(  km) 

e 

w(  *  ) 

hp( km ) 

0 

20000.00000 

.81500000 

35.00000 

302.80000 

120 

20000.32918 

.81479121 

96.77242 

307.03684 

240 

20000.53455 

.81782455 

159.37721 

246.40634 

360 

20000.57169 

.31912787 

225.67056 

220.34600 

480 

20000.42574 

.81753433 

292.68151 

252.19106 

600 

20000.29002 

.81569181 

357.09877 

289.01717 

720 

20000.08480 

.81477000 

60.15629 

307.41565 

Case  14  Continued  (Static  Data) 

_ 

0 

20000.00000 

.81500000 

35.00000 

302.80000 

120 

19998.77302 

.81477669 

96 . 77658 

307.03890 

240 

19994.79723 

.81777109 

159.39778 

246.43019 

J6Q 

19978.02947 

.81892391 

225 . 76558 

220.34343 

480 

19S64. 35529 

.81721454 

292.94967 

251.99383 

600 

19959.31764 

.81533099 

357.59101 

288.66735 

720 

19957 .00771 

.81437821 

60.88218 

mu 

Table  A. 13:  Case  15  Data 


Dynamic  Data 


5000.00000  .27060000 


v(  ‘  ) 


35.00000 


hp( km) 


249.80000 


120 

4999 .55755 

. 26966883 

131.50579 

254.13272 

240 

4999.17581 

.27733025 

226.22172 

215.55313 

360 

4996.54172 

.27670514 

318.84402 

216.77295 

480 

4996.23941 

.26894529 

53.25147 

255.32435 

600 

4996.18743 

.27044997 

149.18642 

247.76869 

720 

4996.09627 

,27790313 

242.90518 

210.46549 

Case  15 

Continued  (Static  Data) 

0 

5000.00000 

.27060000 

35.00000 

249.80000 

120 

4992.72282 

.26863994 

130.99989 

254.27807 

240 

4986.82938 

.27542920 

224.50469 

216.11097 

360 

4981.82717 

.27476373 

315.32094 

215.80178 

480 

4975.91550 

.26627296 

47.01537 

253.76373 

600 

4968.69449 

.26566655 

139.58943 

251.47856 

720 


4962.87185 


27232759  229.45923  214.14492 


Table  A. 14:  Case  16  Data 


I 

I 


Dynamic 

Data 

t 

a(  km ) 

e 

»(  *  ) 

hp( km ) 

0 

10000.00000 

.63530000 

35.00000 

249.80000 

120 

9999.95380 

. 63788457 

338.16373 

223.93753 

24  0 

9999.87033 

.63925371 

281.67942 

210.21608 

360 

9998.83121 

.63836116 

225.28241 

218.76572 

480 

9998.87389 

.63625965 

168.61203 

239.79389 

600 

9998.86252 

.63437761 

11.56627 

258.60798 

720 

9999.01324 

.63428386 

54.41548 

259 . 60056 

Case  16 

Continued  (Static  Data) 

0 

10000.00000 

.63530000 

35.00000 

249.80000 

120 

9991.71527 

.63758789 

338.12500 

223.91864 

240 

9980.55809 

.63856010 

281.48910 

210.17196 

360 

9969.78988 

.63729995 

224.82070 

218.84329 

480 

9959.72721 

.63479035 

167.76473 

240.18852 

600 

9954.98099 

.63274012 

110.24516 

258 . 86516 

720 

9952.39850 

.63264429 

52 . 58695 

258.87037 

Table  A. 

15:  Case  17 

Data 

Dynamic  Data 

nm 

a(  km) 

e 

*C) 

hp( km) 

0 

20000.00000 

.81500000 

35.00000 

302.80000 

120 

20000.06508 

.81675934 

19.03105 

267.62517 

240 

20000.01976 

.81658411 

3.07461 

271.12147 

100 


360 

20000.01611 

.81758658 

347.07064 

251.07133 

•180 

20000 . 12724 

.81638078 

331.20081 

275 . 20775 

600 

20000. 12863 

.81594742 

315.40398 

283.87521 

720 

20000.11430 

.81589036 

299.70516 

285.01381 

Case  l 7 

Continued  (Static  Datal 

0 

20000.00000 

.81500000 

35.00000 

302.80000 

120 

19997.34369 

.81673420 

19.02958 

267.62926 

2-10 

19991.65376 

.81650769 

3.06544 

271.11480 

360 

19985.81237 

.81745654 

34  7 . 04657 

251.07940 

•180 

19979.58406 

.81619343 

331 . 15241 

275 . 17885 

600 

19976.79774 

.81573750 

315.32567 

283.77468 

720 

19974.07560 

.81565400 

299.59326 

284.94085 
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Appendix  B 


Supplemental  Figures  lor  Sections  1  and  ill 

The  following  pages  contain  two  dimensional  graphs  of 
the  data  tabulated  in  Appendix  A. 
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Fig.  B.6:  Semimajor  Axis  vs.  Time  (Case  4) 


'*111 


o 

o 


X}|3IJ}U9303 


109 


120 


o 


IP 

4/ 

O 


<1> 

C 


V) 

> 

“O 


111 


Fig.  B.9:  Periapsis  Alt 


3660- 


<o<0<0intnu*>min 

fOfOfOrOiOfOiOiO 


(uu>j)  sixy  jofDLUiuuas 


112 


Fig.  B.10:  Semimajor  Axis  vs.  Time  (Case 
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Fig.  B.  11:  Eccentricity  vs.  Time  (Case  5) 
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B.12:  Periapsis  Altitude  vs.  Time  (Cci 
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Fig.  B.13:  Semimajor  Axis  vs.  Time  (Case  6) 
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Fig.  B.15:  Periapsis  Altitude  vs.  Time  (Cci 
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Fig.  B.16:  Semimajor  Axis  vs.  Time  (Case  7) 
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B.17:  Eccentricity  vs.  Time  (Case  7) 
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Fig.  B.18:  Periapsis  Altitude  vs.  Time  (Ccis 
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Fig.  B.19:  Semimajor 
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B.20:  Eccentricity  vs.  Time  (Case  S) 
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Fig.  B.21:  Periapsis  Altitude  vs.  Time 
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« 

Sernimajor  Axis  vs.  Time  (Case  9) 
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Fig.  B.24:  Periaps 


Fig.  B.25:  Semimajor  Axis  vs.  Time  (Cas 
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Fig.  B .50:  Periapsis  Altitude  vs.  Time  l1- 
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B.31:  Semimajor  Aas  vs.  Time  (Case 
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Fig.  B.34:  Semimajor 
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Eccentricity  vs.  Time  (Case 


138 


Periaps 


Semimajor  Axis  vs.  Time  (Ca 


|  I  I  I  I  I  I  I  I  I  I  I  I  II  I  I  I  I  I  I  I  I  M  I  I  I  I  II  I  I  I  II  I  I  !  I  I  I  I  I  I  1  I  ■  I  ■  I  ■  I  I  ■  III  I  I  ■  I  M  I  M 


§ 

8 

o 

o 

8 

8 

o 

o 

8 

o 

o 

a 

00 

O) 

CO 

<0 

m 

*fr 

(O 

00 

f** 

00 

00 

5 

00 

5 

5 

Ll_ 

6 

d 

o 

o 

o 

o 

o 

O 

X}j3U}iis3o;3 


140 


•J  J 


w} 

o 

l  > 


<d 

<— 

c 

i — 


(n 

> 

v 

V 

Z3 


10 

'</) 

Cl 

O 

i 

4> 

CL 


cm 

K) 

CD 


g> 

Li- 


141 


142 


Fig.  B.40:  Semimajor 
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Fig.  B.43:  Semimajor 
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Fig.  B.44:  Eccentricity  vs.  Time  (Case  1t*i 
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Fig.  B.45:  Periapsis  Altitude  vs. 


.46:  Semimajor  Axis  vs.  Time  (Case 
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Fig.  B.47:  Eccentricity  vs.  Time  (Co 


4) 

■o 


(lu>i)  gpn^mv  SjsdDusj 


151 


Fig.  B.49:  Periapsis  Altitude  vs.  Time  (ChecI  Casej 
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Fig.  B.5I  Eccentricity  vs.  Time  (Chect  Case) 
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Fig. B. 52:  Atmospheric  Density  vs.  Time 
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